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INTRODUCTION 

The seminal work by Bardeen, Cooper and Schrieffer (BCS) [l] extended further by Eliash- 
berg [3] to the intermediate coupling regime solved one of the major scientific problems of 
Condensed Matter Physics in the last century. While the BCS theory provides a qualita- 
tively correct description of some novel superconductors like magnesium diboride and doped 
fullerenes, if the phonon dressing of carriers (i.e. polaron formation) is properly taken into 
account, high-temperature superconductivity (HTS) of cuprates represents a challenge to 
the conventional theory. Here I discuss a multi-polaron approach to the problem based on 
our extension of the BCS theory to the strong-coupling regime {3] . Attractive electron corre- 
lations, prerequisite to any HTS, are caused by an almost unretarded electron-phonon (e-ph) 
interaction sufficient to overcome the direct Coulomb repulsion in this regime. Low energy 
physics is that of small polarons and bipolarons (real-space electron (hole) pairs dressed by 
phonons). They are itinerant quasiparticles existing in the Bloch states at temperatures 
below the characteristic phonon frequency. Since there is almost no retardation (i.e. no 
Tolmachev-Morel-Anderson logarithm) reducing the Coulomb repulsion, e-ph interactions 
should be relatively strong to overcome the direct Coulomb repulsion, so carriers must be 
polaronic to form pairs in novel superconductors. I identify the Frohlich electron-phonon in- 
teraction as the most essential for pairing in superconducting cuprates. Many experimental 
observations have been satisfactorily understood in the framework of the bipolaron theory 
(3I providing evidence for a novel state of electronic matter in layered cuprates, which is a 
charged Bose-liquid of small mobile bipolarons. 

Here the band structure and essential interactions in oxide superconductors are discussed 



in section 1, and the "Frohlich-Coulomb" model of HTS is introduced in section 2, including 
discussions of single-polaron (2.1, 2.2, 2.3, 2.4, 2.5) and multipolaron (2.4, 2.6, 3.1, 3.2) 
problems, low-energy structures (3.3), and the phase diagram of cuprates (3.4). "Individ- 
ual" versus Cooper pairing (3.5), normal state properties (section 4), in particular in-plane 
resistivity, the Hall effect, magnetic susceptibility and the Lorenz number (4.1), the Nernst 
effect (4.2), diamagnetism (4.3), spin and charge pseudogaps, and c-axis transport (4.4) 
are also discussed. I present a parameter-free evaluation of T c (5.1), and an explanation 
of isotope effects (5.2), specific heat anomaly (5.3), upper critical fields (5.4), symmetries 
and space modulations of the order parameter (5.5), and a model of overdoped cuprates as 
mixtures of mobile bipolarons and degenerate lattice polarons (section 6). 



I. BAND STRUCTURE AND ESSENTIAL INTERACTIONS IN CUPRATES 

A significant fraction of theoretical research in the field of HTS has suggested that the 
interaction in novel superconductors is essentially repulsive and unretarded, and it could 
provide high T c without phonons. Indeed strong on-site repulsive correlations (Hubbard U) 
are essential in shaping the insulating state of undoped (parent) compounds. Different from 
conventional band-structure insulators with completely filled and empty Bloch bands, the 



filled band clS db result of the Coulomb 



4] 



Mott insulator arises from a potentially metallic hal 
blockade of electron tunnelling to neighboring sites 

In our approach to cuprate superconductors we take the view that cuprates and related 
transition metal oxides are charge-transfer Mott-Hubbard insulators at any relevant level of 
doping P|. The one-particle density-of-states (DOS) of cuprates is schematically represented 
by Fig.l, as it has been established in a number of site-selective experiments [5| and in 
the first-principle numerical ("LDA+U") [6] and semi-analytical cluster 7] band structure 
calculations properly taking into account the strong on-site repulsion. Here d-band of the 
transition metal (Cu) is split into the lower and upper Hubbard bands by the on-site repulsive 
interaction U, while the first band to be doped is an oxygen band within the Hubbard gap. 
The oxygen band is less correlated and completely filled in parent insulators, so a single 
oxygen hole has well defined quasi-particle properties in the absence of interactions with 
phonons and with spin fluctuations of d-band electrons. 

Unfortunately, the Hubbard U model shares an inherent difficulty in determining the 
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FIG. 1: DOS in cuprates. The chemical potential \x is inside the charge transfer gap as observed 



in the tunnelling experiments 



221 ] because of bipolaron formation |3( ■ It could enter the oxygen 



band in overdoped cuprates, if bipolarons coexist with unpaired degenerate polarons (section 6). 



order when the Mott-Hubbard insulator is doped. While some groups have claimed that 
it describes high-T c superconductivity at finite doping, other authors could not find any 
superconducting instability. Therefore it has been concluded that models of this kind are 
highly conflicting and confuse the issue by exaggerating the magnetism rather than clarifying 
it |a]. The Hubbard-f/ model of HTS and its strong-coupling "t — J" approximation |2| 
are also refutable on experimental ground. A characteristic magnetic interaction, which is 
allegedly responsible for pairing in the model, is the spin-exchange interaction, J = 4t 2 /U, of 
the order of 0.1 eV (here t is the hopping integral). On the other hand, a simple parameter- 
free estimate of the Frohlich electron-phonon interaction (routinely neglected within the 

n 

Hubbard U approach) yields the effective attraction as high as 1 eV |3j. This estimate is 

obtained us,„g the foliar expression for the polaron level shift, E p , the h lg h-f r e q nency^, 

and the static, e , dielectric constants of the host insulator, measured experimentally [10J, 
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where k -1 = — Cq 1 and the size of the integration region is the Brillouin zone (BZ). Since 
= 5 and eo = 30 in La2Cu04 one obtains E p = 0.65 eV. Hence the attraction, which is 
about 2E P , induced by the long-range lattice deformation in parent cuprates is one order of 
magnitude larger than the exchange magnetic interaction. There is virtually no screening of 
e-ph interactions with c— axis polarized optical phonons in doped cuprates because the upper 
limit for an out-of-plane plasmon frequency (< 200 cm -1 ) 11| is well below characteristic 



phonon frequencies, u ~ 400 - 1000 cm _1 . Hence the Frohlich interaction remains the most 
essential pairing interaction at any doping. 

Further compelling evidence for the strong c-ph interaction has come from isotope ef- 

I I 

'ccts [12], more recent high resolution angle resolved photoemission spectroscopics (ARPES) 
and a number of earlier optical 1J, 15, 16, [l7| and neutron-scattering |l8| studies 



of cuprates. The strong coupling with optical phonons, unambiguously established in all 
high-temperature superconductors, transforms holes into lattice mobile polarons and mobile 



superconducting bipolarons as has been proposed jl9( prior the discovery |20|, |21|. 

When the e-ph interaction binds holes into intersite oxygen bipolarons |3j, the chemical 
potential remains pinned inside the charge transfer gap. It is found at a half of the bipolaron 
binding energy, Fig.l, above the oxygen band edge shifted by the polaron level shift E p , 
as clearly observed in the tunnelling experiments by Bozovic et al. in optimally doped 



Lax 85 Sr 15 Cu 4 |22[. The bipolaron binding energy as well as the singlet-triplet bipolaron 
exchange energy (section 3) are thought to be the orig in of normal state charge and spin 
pseudogaps, respectively, as has been proposed by us |23J and later found experimentally [24] . 
In overdoped samples carriers screen part of the e-ph interaction with low frequency phonons. 
Hence, the bipolaron binding energy decreases |25j and the hole bandwidth increases with 
doping. As a result, the chemical potential could enter the oxygen band in overdoped samples 
because of an overlap of the bipolaron and polaron bands, so a Fermi-level crossing could 
be seen in ARPES (section 6). 

II. "FROHLICH-COULOMB" MODEL OF HTS 



Experimental facts tell us that any realistic description of high temperature supercon- 
ductivity should treat the long-range Coulomb and unscreened e-ph interactions on an e qual 
footing. In the past decade we have developed a "Frohlich-Coulomb" model (FCM) 0,13,13 
to deal with the strong long-range Coulomb and the strong long-range e-ph interactions in 
cuprates and other related compounds. The model Hamiltonian explicitly includes a long- 
range electron-phonon and the Coulomb interactions as well as the kinetic and deformation 
energies. The implicitly present large Hubbard U term prohibits double occupancy and 
removes the need to distinguish fermionic spins since the exchange interaction is negligible 
compared with the direct Coulomb and the electron-phonon interactions. 



Introducing spinless fermionic, c n , and phononic, d ma , operators the Hamiltonian of the 
model is written as 



H = - ^2 [ T ( n ~ n ') c n c n' - V c (n - n')4c n c^,c n / 
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ma 

where T(n) is the hopping integral in a rigid lattice, e a is the polarization vector of the 
ath vibration coordinate, u m _ n = (m — n)/|m — n| is the unit vector in the direction from 
electron n to ion m, g a (m — n) is the dimensionless e-ph coupling function, and V c (n — n') 
is the inter-site Coulomb repulsion. g a (m — n) is proportional to the force acting between 
the electron on site n and the ion on m. For simplicity, we assume that all the phonon 
modes are non- dispersive with the frequency u a . We also use h = k B = c = 1. 

If the electron-phonon interaction is strong, i.e. the conventional e-ph coupling constant 
of the BCS theory is large, A„ >4j then the weak-conp.ing BCS Q and the intermediate- 



coupling Migdal-Eliashberg [2|, |28( approaches cannot be applied [29]. Nevertheless the 



Hamiltonian, Eq.(2), can be solved analytically by using the "1/A" multi-polaron expansion 
technique 0, if A = E p /zT(a) > 1. Here the polaron level shift is 

E p = ^2uj a g 2 a (n)(e a -u n ) 2 , (3) 
lid 

and zT(a) is about the half-bandwidth in a rigid lattice. As I discuss below, the model 
shows a rich phase diagram de pen ding on the ratio of the inter-site Coulomb repulsion V c 
and the polaron level shift E p |27|. The ground state of FCM is a polaronic Fermi liquid 
when the Coulomb repulsion is large, a bipolaronic high-temperature superconductor at in- 
termediate Coulomb repulsions, and a charge-segregated insulator if the repulsion is weak. 
FCM predicts superlight polarons and bipolarons in cuprates with a remarkably high su- 
perconducting critical temperature. Cuprate bipolarons are relatively light because they are 
inter-site rather than on-site pairs due to the strong on-site repulsion, and because mainly 
c-axis polarized optical phonons are responsible for the in-plane mass renormalization. The 
relatively small mass renormalization of polaronic and bipolaronic carries in FCM has been 
confirmed numerically using the exact QMC [3^ , cluster diagonalization 
32j simulations. 
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(Bi)polarons describe many properties of cuprates [3j], in particular normal-state trans- 
port (section 4), including in-plane and out-of-plane resistivity, the Hall effect, spin sus- 
ceptibility, thermal conductivity, normal state pseudogaps, the Nernst effect, normal state 
diamagnetism, superconducting transition, including high values of T c , isotope effects, un- 
usual upper critical fields, symmetries and real-space modulations of the superconducting 
order parameter (section 5). 



A. Single lattice polaron 

Let us first discuss a single lattice-polaron problem. Conducting electrons in inorganic 
and organic matter interact with vibrating ions. If phonon frequencies are sufficiently low, 
the local deformation of ions, caused by electron itself, creates a potential well, which traps 
the electron even in a perfect crystal lattice. This self-tray ping phenomenon was predicted 
Landau 0L ^ was studied in greater detail by Pekar J^ij], Frohlich j^, Feynman 
"\, Rashba \%*\. Devreesej^] and other authors in the effective mass approximation for 
the electron placed in a continuous polarizable medium, which leads to a so-called large 
or continuous polaron. Large polaron wave functions and corresponding lattice distortions 
spread over many lattice sites. The trapping is never complete in the perfect lattice. Due to 
finite phonon frequencies ion polarizations follow polaron motion if the motion is sufficiently 
slow. Hence, large polarons with a low kinetic energy propagate through the lattice as free 
electrons but with an enhanced effective mass. 

When the electron-phonon (e-ph) interaction energy E p is compared with the electron 
energy-bandwidth, all electrons in the Bloch bands of the crystal are "dressed" by phonons. 
In this strong-coupling regime, A = E p /D > I, the finite bandwidth 2D becomes important, 
so the continuous approximation cannot be applied. The main features of small polarons 
were understood by Tjablikov jsil, Y amashita and Kurosava (^J , Sewell j^lj, Holstein j3| 
and his school, Lang and Firsov Eagles [d^ , and by other researches and described 
in several review papers and textbooks j45|. The polaron shift of the atomic level and an 
exponential reduction of the bandwidth (see below) at large values of A are among those 
features. The shift can be easily understood using a toy model of an electron localized on 
site n and interacting with a single ion vibrating near site m in the direction connecting n 



and m, Fig.2. The vibration part of the Hamiltonian in this toy model is 

Id 2 kx 2 

where M is the ion mass, k = Muj 2 is the spring constant, and x is the ion displacement. 
The electron potential energy due to its Coulomb interaction with the ion is approximately 

V = V Q (l-x/a), (5) 

where Vq = —Ze 2 /a is the Coulomb energy in a rigid lattice (an analog of the crystal field 
potential), and a is the average distance between sites. Hence the Hamiltonian of the model 
is given by 

H = E^ + fM-—^ + —, (6) 

where E a is the atomic level at site m in the rigid lattice, which includes the crystal field, 
/ = Ze 2 /a 2 is the Coulomb force, and n = c^c is the occupation number operator on site n 
expressed in terms of the electron annihilation c and creation c* operators. This Hamiltonian 
can be readily diagonals using a displacement transformation of the vibration coordinate x, 

x = y-hf/k. (7) 
The transformed Hamiltonian has no electron-phonon coupling, 



where we used h 2 = h because of the Fermi statistics. It describes a small polaron at the 
atomic level, shifted by the polaron level shift E p = f 2 /2k, and entirely decoupled from 
ion vibrations. The ion vibrates near a new equilibrium, shifted by f/k, with the "old" 
frequency u. As a result of the local ion deformation, the total energy of the whole system 
decreases by E p since a decrease of the electron energy by —2E P overruns an increase of the 
deformation energy E p . 

The tunnelling of small polarons in the lattice can be understood within a simple Holstein 



model 42j consisting of two molecules and a single electron. Here I slightly simplify the 
original Holstein model replacing two molecules by two rigid sites 1 ("left") and 2 ("right") 
with the hopping amplitude t between them. The electron interacts with a vibrational mode 




FIG. 2: Localized electron shifts the equilibrium position of an ion and lowers the atomic energy 
level. 

of an ion, placed at some distance in between, Fig. 3, rather than with the intra-molecular 
vibrations: 

H = t(c\c 2 + eld) + H ph + H e _ ph , (9) 

where H e _ ph depends on the polarization of vibrations, and E a = is taken. If the ion 
vibrates along the perpendicular direction to the hopping (in "c" -direction) we have 

H e ^p h = fcx(c\cx + 4c 2 ), (10) 

and 

He-ph = faX{c{ci - 4c 2 ), (11) 

if the ion vibrates along the hopping ("a" direction). 

The wave-function of the electron and the ion is a linear superposition of two terms 
describing the electron on the "left" and on the "right" site, respectively, 

^ = [u(x)c\ + v(x)4}\0), (12) 

where |0) is the vacuum state describing a rigid lattice without the extra electron. Substi- 
tuting ip into the Schrodinger equation, Hip = Eip, we obtain two coupled equations for the 



-< — ► a-polarised 
c-polarisedj (^^^j 

h o 

"left" "right" 



FIG. 3: Electron tunnels between sites 1 ("left") and 2 ("right") with the amplitude t and interacts 
with c-axis or a-axis polarized vibrational modes of the ion, placed in between. 

amplitudes, 

{E- f a , c x-H ph )u{x) = tv{x), (13) 

(E ± f a , c x - H ph )v(x) = tu(x). (14) 

There is the exact solution for the "c"-axis polarization, when a change in the ion position 
leads to the same shift of the electron energy on the left and on the right sites, 

u(x) = uxn(x), 

V (x) = VXn(x), (15) 

where u and v are constants and 

Xn(x) = (^^y\^ (16) 

is the harmonic oscillator wave-function. There are two ladders of levels given by 

E± = -E p ± t + oo(n + 1/2) (17) 



with E p = pj2k. Here 

H n (0 = (-i) B ^^r (is) 

are the Hermite polynomials, and n — 0, 1, 2, 3, ... . Hence the c-axis single-ion deformation 
leads to the polaron level shift but without any renormalization of the hopping integral 
t. In contrast, a-polarized vibrations with the opposite shift of the electron energy on the 
left and on the right sites, strongly renormalize the hopping integral. There is no simple 
general solution of the Holstein model in this case, but one can find it in two limiting cases, 
nonadiabatic, when t <^ u and adiabatic, when t ^> uj. 



B. Non-adiabatic small polaron 



In the non-adiabatic regime the ion vibrations are fast and the electron hopping is slow. 
Hence one can apply a perturbation theory in powers of t to solve 



E-faX- H ph -t 

-t E + f a x - H ph 




0. 



(19) 



We take t = in zero order, and obtain a two-fold degenerate ground state [u l ' r (x),v l,r (x)], 
corresponding to the polaron localized on the left (/) or on the right (r) sites, 

Mu 



u l (x) 



exp 



(x + fa/ky 



v\x) = 



(20) 



and 



v r (x) 



exp 



Mu 



u r {x) = 0, 

X - fa/kf 



(21) 



with the energy Eq = —E p + u/2, where E p = f% /2k. The eigenstates are found as linear 
superpositions of two unperturbed states, 
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v r (x) 



(22) 



Here the coefficients a and (5 are independent of x. The conventional secular equation for E 
is obtained, multiplying the first row by u\x) and the second row by v r (x), and integrating 



over the vibration coordinate, x, each of two equations of the system. The result is 

f E-E t \ 
det =0 (23) 

^ t E-E J 

with the renormalized hopping integral 

t SZo dx \ ul ( x )\ 2 

The corresponding eigenvalues, E± are 

E± = uj/2 -E p ± t. (25) 

The hopping integral splits the degenerate level, as in the rigid lattice, but an effective 
'bandwidth' 2t is significantly reduced compared with the bare one 

t = texp(-2E p /Lu). (26) 

This polaron band narrowing originates in a small overlap integral of two displaced oscillator 
wave functions u l (x) and v r (x). 

C. Adiabatic small polaron 

In the adiabatic regime, when t ^> u, the electron tunnelling is fast compared with the 
ion motion. Hence one can apply the Born-Oppenheimer adiabatic approximation taking 
the wave function in the form 



(27) 



Here u a (x) and v a (x) are the electron wave functions obeying the Schrodinger equation with 
the frozen ion deformation x, i.e. 
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E a {x) ~ f a X ~t \ ( U a (x) 

-t E a (x) + f a x \ v a (x) 



= 0. (28) 



The lowest energy level is found as 

E a (x) = - v /(f a xy + ti. (29) 



E a (x) together with kx 2 /2 play the role of a potential energy term in the equation for the 
'vibration' wave function, x( x )i 



1 d 2 kx 2 



+ ^r-V(fax) 2 + t 2 



X (x) = E X (x). (30) 



2Mdx 2 2 

Terms with the first and second derivatives of the electron wave-functions u a (x) and v a (x) are 
small compared with the corresponding derivatives of X ( x ) m the adiabatic approximation, so 
they can be neglected in Eq.(31). As a result we arrive with the familiar double-well potential 



problem, where the potential energy U(x) = kx 2 /2 — a/ (fa x ) 2 + 1 2 has two symmetric 
minima, separated by a barrier. Minima are located approximately at 

x m = ±f a /k (31) 

in the strong-coupling limit, E p ^> t, and the potential energy near the bottom of each 
potential well is about 

U(x) = -E p + K\A-Ulk)\ (32) 

If the barrier were impenetrable, there would be the ground state energy level E = —E p + 
uj/2, the same for both wells. The tunnelling under the barrier results in a splitting of this 
level 2t, which corresponds to a polaron band in the lattice. It can be estimated using the 
quasi-classical approximation as 



t oc exp 



2 / p(x)dx 
Jo 



(33) 



where p(x) = ^2M[U(x) — E ] w (Mk) l l 2 \x — f a /k\ is the classical momentum 
Calculating the integral one finds the exponential reduction of the "bandwidth", 

t oc exp(-2E p /iu), (34) 

which is the same as in the nonadiabatic regime. Holstein found also corrections to this ex- 
pression up to terms of the order of 1/A 2 , which allowed him to estimate the pre-exponential 
factor as 



t ss y/E p ujexp(-2E p /uj). (35) 

The term in front of the exponent differs from t of the non-adiabatic case. It is thus 
apparent that the perturbation approach covers only part of the entire lattice polaron region, 
A > 1. The upper limit of applicability of the perturbation theory is given by t < sJE p uj. 
For the remainder of the region the adiabatic approximation is more accurate. 



D. "1/A" expansion technique: polaron band 



The kinetic energy is small compared with the interaction energy as long as A > 1. Hence 
an analytical approach to the multi-polaron problem is possible with the "1/A" expansion 
technique jjjj], which treats the kinetic energy as a perturbation. The technique is based on 
the fact, known for a long time, that there is an analytical exact solution of a single polaron 
problem in the strong-coupling limit A — > oo. Following Lang and Firsov j^J we apply the 
canonical transformation e s diagonalizing the Hamiltonian, Eq.(2). The diagonalization is 
exact, if T(m) = (or A = oo). In the Wannier representation for electrons and phonons, 

S = g a (m - n)(e Q • u m _ n )4c n (d ma - d ma ). 

m^n,a 

The transformed Hamiltonian is 

H = e- s He s = ^nn'C^c n , + u ^ (d ] ma d ma + - J + (36) 

n^n' ma ^ ' 

} J v(n - n')4c n c^c n / - E p ^ c|,c n , 

n^n' n 

where for simplicity we take u a = u. The last term describes the energy gained by polarons 
due to the e-ph interaction. The third term on the right-hand side is the polaron-polaron 
interaction, 

v(n - n') = V c (n - n') - V ph (n - n'), (37) 

where 

V ph (n - n) = 2u g a (m - n)g a (m. - n) x 

m,a 

(gq, • u m _ n )(eQ, • u m _ n '). 

The phonon-induced interaction V p h is due to displacements of common ions caused by two 
electrons. Finally, the transformed hopping operator cr nn / is given by 



Onn' = ^( n _ n ') exp 



y^[flr a (m-n) (e Q • u m _ n ) (38) 



- g a (m- n')(e Q ■ u m _ n /)] (S ma - 

This term is perturbation at a large A. It accounts for the polaron and bipolaron tunnelling 
and high temperature superconductivity [jj . In particular crystal structures like perovskites, 



a bipolaron tunnelling appears already in the first order in T(n) (see below), so that <j nn / 
can be averaged over phonon vacuum, if the temperature is low enough, T uj. The result 
is 



t(n - n') = ((a nn ,)) vh = T(n - n') exp[-<r> - n')], (39) 



where 



g 2 (n-n') = g a (m - n) (e a • u m _ n ) 



[g a (m - n)(e Q • u m _ n ) -g a (m- n')(e Q ■ u m _ n /)] . 

By comparing Eqs.(40) and Eqs.(38,3), the bandwidth renormalization exponent can be 
expressed via E p and V p h as follows 



g 2 (n-n') = - 

10 



E p -±V ph (n-n') 



(40) 



In zero order with respect to the hopping the Hamiltonian, Eq.(37) describes localized 
polarons and independent phonons, which are vibrations of ions around new equilibrium 
positions depending on the polaron occupation numbers. The phonon frequencies remain 
unchanged in this limit. The middle of the electron band falls by the polaron level-shift E p 
due to a potential well created by lattice deformation. The finite hopping term leads to the 
polaron tunnelling because of degeneracy of the zero order Hamiltonian with respect to site 
positions of the polaron. 



E. From continuous to small Holstein and small Frohlich polarons: QMC simula- 
tion 

The narrowing of the band and the polaron effective mass strongly depend on the ra- 



dius of the electron-phonon interaction [26[. Let us compare the small Holstein polaron 
(SHP) formed by a short-range e-ph interaction and a small polaron formed by a long-range 
(Frohlich) interaction, which we refer to as the small Frohlich polaron (SFP). For simplicity 
we consider the interaction with a single phonon branch. In general, there is no simple rela- 
tion between the polaron level-shift E p and the exponent g 2 . This relation depends on the 
form of the electron-phonon interaction. In the nearest-neighbor approximation the effective 
mass renormalization is given by 

m* /m = e ff2 , 
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FIG. 4: One-dimensional model of the lattice polaron on the chain interacting with displacements 
of all ions of another chain. 

where m is the bare band mass and g 2 = g 2 (a). 

If the interaction is short-ranged, g a (m) = K5 m p (the Holstein model), then g 2 = E p ju. 



than 1. To estimate 7 let us consider a one- dimensional chain model with the long-range 
Coulomb interaction between the electron on chain ( x ) and ion vibrations of another chain 
(o), polarized in the direction perpendicular to the chains (3c]], Fig.4. The corresponding 
force is given by 



Here the distance along the chains |m — n| is measured in units of the lattice constant a, the 
inter-chain distance is also a, and we take a = 1. For this long-range interaction we obtain 
7 = 0.39. Hence the effective mass renormalization is much smaller than in the Holstein 



Not only does the small polaron mass strongly depend on the radius of the electron- 
phonon interaction, but also does the range of the applicability of the analytical 1/ A expan- 
sion theory. The theory appears almost exact in a wide region of parameters for the Frohlich 
interaction. The exact polaron mass in a wide region of the adiabatic parameter u/T(a) 
and coupling was calculated with the continuous-time path-integral Quantum Monte Carlo 
(QMC) algorithm Q. This method is free from any systematic finite-size, finite-time-step 
and finite-temperature errors and allows for an exact (in the QMC sense) calculation of the 



Here k is a constant. In general, we have g 2 



^E p /uo with the numerical coefficient 7 less 




(41) 





FIG. 5: Inverse effective polaron mass in units of 1/m = 2T(a)a 2 (T(a) = t on the graph) |3(| 



ground-state energy and the effective mass of the lattice polaron for any electron-phonon 
interaction. At large A (> 1.5) SFP was found to be much lighter than SHP, while the large 
Frohlich polaron (i.e. at A < 1) was heavier than the large Holstein polaron with the same 
binding energy, Fig. 5. The mass ratio m* FP /m* HP is a non-monotonic function of A. The 
effective mass of the Frohlich polaron, m FP (X) is well fitted by a single exponent, which 
is e 073A for uo = T(a) and e 1,4A for uo = 0.5 T(a). The exponents are remarkably close to 
those obtained with the Lang-Firsov transformation, e°' 78A and e L56A , respectively. Hence, 
in the case of the Frohlich interaction the transformation is perfectly accurate even in the 
moderate adiabatic regime, u/T(a) < 1 for any coupling strength. It is not the case for the 
Holstein polaron. If the interaction is short-ranged, the same analytical technique is applied 
only in the nonadiabatic regime u/T(a) > 1. 

Another interesting point is that the size of SFP and the length, over which the distortion 
spreads, are dif f event. In the strong-coupling limit the polaron is almost localized on one 
site m. Hence, the size of its wave function is the atomic size. On the other hand, the 
ion displacements, proportional to the displacement force g a (m — n), spread over a large 



distance. Their amplitude at a site n falls with the distance as |m— n|~ 3 in our one- 
dimensional model. The polaron cloud (i.e. lattice distortion) is more extended than the 
polaron itself. Such polaron tunnels with a larger probability than the Holstein polaron due 
to a smaller relative lattice distortion around two neighboring sites. For the short-range 
e-ph interaction the entire lattice deformation disappears at one site and then forms at its 
neighbor, when the polaron tunnels from site to site. Therefore 7 = 1 and the polaron is 
very heavy already at A ~ 1. On the contrary, if the interaction is long-ranged, only a 
fraction of the total deformation changes every time the polaron tunnels from one site to its 
neighbor, and 7 is smaller than 1. A lighter mass of SFP compared with the nondispersive 
SHP is a generic feature of any dispersive electron-phonon interaction. 



F. Attractive correlations of small polarons 



Lattice deformation also strongly affects the interaction between electrons. At large 
distances polarons repel each other in ionic crystals, but their Coulomb repulsion is sub- 
stantially reduced due to ion polarization. Nevertheless two large polarons can be bound 
into a large bipolaron by an exchange interaction even with no additional e-ph interaction 



but the Frohlich one 
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When a short-range deformation potential and molecular e-ph interactions (e.g. of the 
Jahn- Teller type are taken into account together with the long-range Frohlich inter- 
action, they can overcome the Coulomb repulsion j^J. The resulting interaction becomes 
attractive at a short distance of about a lattice constant. Then two small polarons readily 
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54, 



55j | . because their band is nar- 



form a bound state, i.e. a small bipolaron 
row. Consideration of particular lattice structures shows that small bipolarons are mobile 
even when the electron-phonon coupling is strong and the bipolaron binding energy is large 
(3] (see below). Here we encounter a novel electronic state of matter, a charged Bose liquid 
of electron molecules with double elementary charge 2e, qualitatively different from normal 
Fermi-liquids in ordinary metals and from the Bardeen-Cooper-Schrieffer (BCS) superffuids 
in conventional superconductors. 

The origin of the attractive force between two small polarons can be readily understood 
from about the same toy model as in Fig. 2, but with two electrons on neighbor sites 1,2 
interacting with an ion in between 3, Fig. 6. For generality we now assume that the ion is a 




FIG. 6: Two localized electrons shift the equilibrium position of the ion (3). As a result two 
electrons on neighboring sites 1 and 2 attract each other. 

three-dimensional oscillator described by a displacement vector u, rather than by a single- 
component displacement x as in Fig. 2. The vibration part of the Hamiltonian in the model 
is 

1 / d \ 2 ku 2 
^ = -2mU) (42) 

Electron potential energies due to the Coulomb interaction with the ion are approximately 

V 1 , 2 = V (l-u-e 1 , 2 /a), (43) 

where e 12 are units vectors connecting sites 1, 2 and site 3, respectively. Hence the Hamil- 
tonian of the model is given by 

1 / d \ 2 ku 2 

H = £ a (ni + n 2 ) + u • ftnx + f 2 n 2 ) - — ^— J + — , (44) 

where fi j2 = ^e 2 e 12 /a 2 is the Coulomb force, and n 12 are occupation number operators at 
every site. This Hamiltonian is also readily diagonalized by the same displacement trans- 
formation of the vibronic coordinate u as above, 

u = v - (fxni + f 2 n 2 ) jk. (45) 



The transformed Hamiltonian has no electron-phonon coupling, 

1 / d \ 2 fcf 2 

H=(E a - E p )(h! + n 2 ) + V vh h x h 2 - — I- + — , (46) 

and describes two small polarons at their atomic levels shifted by the polaron level shift 
E p = ff 2 /2k, which are entirely decoupled from ion vibrations. As a result, the lattice 
deformation caused by two electrons leads to an effective interaction between them, V p h, 
which should be added to their Coulomb repulsion, V c , 

V ph = -fx • h/k. (47) 

When V p h is negative and larger by magnitude than the positive V c the resulting interaction 
becomes attractive. 

Applying the polaron canonical transformation to a generic "Frohlich-Coulomb" Hamil- 



tonian, allows us explicitly calculate the effective attraction of small polarons |27(, Eq.(38), 
and elaborate more physics behind the lattice sums in Eq.(3,38) and Eq.(41). If a carrier 
(electron or hole) acts on an ion with a force f, it displaces the ion by some vector x = f/k. 
Here k is the ion's force constant. The total energy of the carrier-ion pair is — f 2 /(2k). 
This is precisely the summand in Eq.(3) expressed via dimensionless coupling constants. 
Now consider two carriers interacting with the same ion, see Fig. 7. The ion displacement 
is x = (Fi + F 2 )//c and the energy is — fj 2 /(2/c) — F|/(2fc) — (fi ■ f 2 )/fc. Here the last term 
should be interpreted as an ion-mediated interaction between the two carriers. It depends 
on the scalar product of fi and f 2 and consequently on the relative positions of the carriers 
with respect to the ion. If the ion is an isotropic harmonic oscillator, as we assume here, 
then the following simple rule applies. If the angle between fx and f 2 is less than 7r/2 the 
polaron-polaron interaction will be attractive, if otherwise it will be repulsive. In general, 
some ions will generate attraction, and some repulsion between polarons, Fig. 7. The overall 
sign and magnitude of the interaction is given by the lattice sum in Eq.(38), the evaluation 
of which is elementary. One should also note that according to Eq.(41) an attractive in- 
teraction reduces the polaron mass (and consequently the bipolaron mass), while repulsive 
interaction enhances the mass. 




FIG. 7: Mechanism of the polaron-polaron interaction, (a) Together, two polarons (solid circles) 
deform the lattice more effectively than separately. An effective attraction occurs when the angle 
<j> is less than ir/2 . (b) A mixed situation. Ion 1 results in repulsion between two polarons while 
ion 2 results in attraction. 

III. SUPERLIGHT BIPOLARONS IN HIGH-T C CUPRATES 

Consideration of particular lattice structures shows that small inter-site bipolarons are 
perfectly mobile even when the electron-phonon coupling is strong and the bipolaron bind- 
ing energy is large. Let us analyze the important case of copper-based high-T c oxides. As 
discussed in the introduction they are doped charged-transfer ionic insulators with narrow 
electron bands. Therefore, the interaction between holes can be analyzed using computer 
simulation techniques based on a minimization of the ground state energy of an ionic in- 
sulator with two holes, the lattice deformations and the Coulomb repulsion fully taken 
into account, but neglecting the kinetic energy terms. Using these techniques net inter- 
site interactions of the in-plane oxygen hole with the apex hole, Fig. 8, and of two in-plane 
oxygen holes, Fig. 10, were found to be attractive in La 2 Cu0 4 with the binding ener- 
gies A = 119meV and A = 60meV, respectively. All other interactions were found to be 
repulsive. 



FIG. 8: Apex bipolaron tunnelling in perovskites 



A. Apex bipolarons 

Both apex and in-plane bipolarons can tunnel from one unit cell to another via the 
smg/e-polaron tunnelling from one apex oxygen to its apex neighbor in case of the apex 
bipolaron j^J, Fig.8, or via the next-neighbor hopping in case of the in-plane bipolaron [27], 
Fig. 10. The Bloch bands of these bipolarons are obtained using the canonical transformation, 
described above, projecting the transformed Hamiltonian, Eq.(37), onto a reduced Hilbert 
space containing only empty or doubly occupied elementary cells, and averaging the result 
with respect to phonons |3[. The wave function of the apex bipolaron localized, let us say, 
in the cell m is written as 

4 

\ m ) = J2Mcl pex \0), (48) 

i=l 

where i denotes the p XjV orbitals and spins of the four plane oxygen ions in the cell, Fig.8 
and c\ pex is the creation operator for the hole on one of the three apex oxygen orbitals with 
the spin, which is the same or opposite of the spin of the in-plane hole depending on the 
total spin of the bipolaron. The probability amplitudes Aj are normalized by the condition 
\Ai\ = 1/2, if four plane orbitals p x i,P y 2,Px3 and p y 4 are involved, or by \Aj\ = l/y/2 if only 
two of them are relevant. Then a matrix element of the Hamiltonian Eq.(37) describing the 
bipolaron tunnelling to the nearest neighbor cell m + a is found as 

t = (m|#|m + a) = ^T^e^ , (49) 

where T^ x e~ 92 is a single polaron hopping integral between two apex ions. The inter- 
site bipolaron tunnelling appears already in the first order with respect to the single-hole 



FIG. 9: Simplified model of the copper-oxygen perovskite layer. 



transfer T^, ex , and the bipolaron energy spectrum consists of two bands E x,y (K), formed 
by the overlap of p x and p y apex oxygen orbitals, respectively: 

E*(K) = tcos(K x ) -t'cos(Ky), (50) 
E y (K) = -t'cos(K x )+tcos(K y ). 

They transform into one another under ir/2 rotation. If t, t' > 0, "x" bipolaron band has its 
minima at K =(±7r, 0) and y-band at K = (0, ±7r). In these equations t is the renormalized 
hopping integral between p orbitals of the same symmetry elongated in the direction of the 
hopping (ppa) and t' is the renormalized hopping integral in the perpendicular direction 
(ppix). Their ratio t/t' = T^ x /T'^f x = 4 as follows from the tables of hopping integrals in 
solids. Two different bands are not mixed because T" p !? = for the nearest neighbors. A 
random potential does not mix them either, if it varies smoothly on the lattice scale. Hence, 
we can distinguish l x' and '?/' bipolarons with a lighter effective mass in x or y direction, 
respectively. The apex z bipolaron, if formed, is ca. four times less mobile than x and y 
bipolarons. The bipolaron bandwidth is of the same order as the polaron one, which is a 
specific feature of inter-site bipolarons. For a large part of the Brillouin zone near (0, tt) for 
V and (tt, 0) for '?/' bipolarons, one can adopt the effective mass approximation 

K 2 K 2 

E x 'y(K) = x — + y — (51) 

9m*' 9<m** 
z, '"'x,y 

with K XtV taken relative to the band bottom positions and m** = 1/t, m** = Am** . 



B. In-plane bipolarons 



Now let us consider in-plane bipolarons in a two-dimensional lattice of ideal octahedra 
that can be regarded as a simplified model of the copper-oxygen perovskite layer, Fig. 9 [27]. 
The lattice period is a = 1 and the distance between the apical sites and the central plane 
is h — a/2 = 0.5. For mathematical transparency we assume that all in-plane atoms, both 
copper and oxygen, are static but apex oxygens are independent three-dimensional isotropic 
harmonic oscillators. Due to poor screening, the hole-apex interaction is purely coulombic, 

g a (m - n) 



m — n 



where a = x,y, z. To account for the experimental fact that z-polarized phonons couple 
to the holes stronger than others we choose /a/2. The direct hole- hole 

repulsion is 

VJn - n'l 1 ' 



V2\n-n'\ 

so that the repulsion between two holes in the nearest neighbor (NN) configuration is V c . 
We also include the bare NN hopping T^n, the next nearest neighbor (NNN) hopping across 
copper Tnnn and the NNN hopping between the pyramids T' NNN . 

The polaron shift is given by the lattice sum Eq.(3), which after summation over polar- 
izations yields 



2 ( 1 h 2 \ , 

E p = 2k x uj q } - + 1 = 31.15k x u , 

p * ^\ m-n m-n 6 / 

m * 1 1 1 1 * 



(52) 



where the factor 2 accounts for two layers of apical sites. For reference, the Cartesian 
coordinates are n = (n x + l/2, n y + 1/2, 0), m = (m x , m y , h), and integers. 
The polaron-polaron attraction is 

m — n'Mm — n| 



,n 9 h 2 + (m — n') • (m — n) . 
V ph n - n') = 4ljk* £ — ; 1 (53) 



Performing the lattice summations for the NN, NNN, and NNN' configurations one finds 
Vph = 1.23 E p , 0.80 E p , and 0.82^, respectively. As a result, we obtain a net inter-polaron 
interaction as v NN — V c — 1.23 E p , v NNN = ^| — 0.80 E p , v' NNN — ^ — 0.82^, and the 
mass renormalization exponents as g 2 NN = 0.38(E p /u), g 2 NNN = 0.60(E p /uj) and (q'nnn) 2 = 
Q.59(E p /uj). 



FIG. 10: Four degenerate in-plane bipolaron configurations A, B, C, and D . Some single-polaron 
hoppings are indicated by arrows. 

Let us now discuss different regimes of the model. At V c > 1.23 E p , no bipolarons are 
formed and the system is a polaronic Fermi liquid. Polarons tunnel in the square lattice 
with t = T NN exp(-0. 38 E p /lu) and t' = T NNN exp(-0.60£ p Av) for NN and NNN hoppings, 
respectively. Since g 2 NNN ~ {q'nnn) 2 one can neglect the difference between NNN hoppings 
within and between the octahedra. A single polaron spectrum is therefore 

E 1 (k) = -E p - 2t'[cosk x + cos k y \ - Atcos(k x /2) cos(k y /2). (54) 

The polaron mass is m* = l/(t + 2t'). Since in general t > t', the mass is mostly determined 
by the NN hopping amplitude t. 

If V c < 1.23 E p then intersite NN bipolarons form. The bipolarons tunnel in the plane 
via four resonating (degenerate) configurations A, B, C, and D, as shown in Fig. 10. In the 
first order of the renormalized hopping integral, one should retain only these lowest energy 
configurations and discard all the processes that involve configurations with higher energies. 
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FIG. 11: FCM phase diagram. 

The result of such a projection is the bipolaron Hamiltonian 

H b = (V e - 3.23 E p ) Y^AiAi + B}Bi + C\d + d\D x ] (55) 

i 

-t' Y^i A l B i + B \ c \ + c \ A + A^i + ^-c] 
i 

n 

where 1 numbers octahedra rather than individual sites, x = (1, 0), and y = (0, 1). A Fourier 
transformation and diagonalization of a 4 x 4 matrix yields the bipolaron spectrum: 



E 2 (K) = V C - 3.23E P ± 2t'[cos(K x /2) ± cos(AT y /2)]. 



(56) 



There are four bipolaronic subbands combined in the band of the width 8t'. The effective 
mass of the lowest band is m** = 2/t'. The bipolaron binding energy is A w 1.23E P — V c . 
Inter-site bipolarons already move in the first order of the single polaron hopping. This 
remarkable property is entirely due to the strong on-site repulsion and long-range electron- 
phonon interactions that leads to a non-trivial connectivity of the lattice. This fact combines 
with a weak renormalization of t' yielding a superlight bipolaron with the mass m** oc 
exp(0.60 E p /u). We recall that in the Holstein model m** cx exp(2E p /uj) jlfll ]. Thus the 
mass of the Frohlich bipolaron in the perovskite layer scales approximately as a cubic root 
of that of the Holstein bipolaron. At even stronger e-ph interaction, V c < 1.1QE P , NNN 



bipolarons become stable. More importantly, holes can now form 3- and 4-particle clusters. 
The dominance of the potential energy over kinetic in the transformed Hamiltonian enables 
us to readily investigate these many-polaron cases. Three holes placed within one oxygen 
square have four degenerate states with the energy 2(V c —1.23E p )+V c / ^/2—0.80E p . The first- 
order polaron hopping processes mix the states resulting in a ground state linear combination 
with the energy E 3 = 2.71V C — 3.2QE p — v 'At 2 + t' 2 . It is essential that between the squares 
such triads could move only in higher orders of polaron hopping. In the first order, they are 
immobile. A cluster of four holes has only one state within a square of oxygen atoms. Its 
energy is E A = 4{V C - 1.23E P ) + 2(^/^-0.80^) = 5.41V; - 6.52£ p . This cluster, as well 
as all bigger ones, are also immobile in the first order of polaron hopping. We would like to 
stress that at distances much larger than the lattice constant the polaron-polaron interaction 
is always repulsive, and the formation of infinite clusters, stripes or strings is prohibited. 
We conclude that at V c < 1.16E P the system quickly becomes a charge segregated insulator, 
Fig.ll. 

The fact that within the window, 1.16E P < V c < 1.23E P , there are no three or more 
polaron bound states, means that bipolarons repel each other. The system is effectively a 
charged Bose-gas, which is a superconductor. The superconductivity window, that we have 
found, is quite narrow. This indicates that the superconducting state in cuprates requires a 
rather fine balance between electronic and ionic interactions, Fig.ll. 



C. Low-energy (bi) polaron energy structure of cuprates 

The considerations set out above lead us to a simple model of cuprates The main 

assumption is that all electrons are bound into small singlet and triplet inter site bipolarons 
stabilized by e-ph interactions. As the undoped plane has a half-filled Cu3d 9 band there 
are no empty states for bipolarons to move if they are inter-site. Their Brillouin zone is 
half of the original electron BZ and is completely filled with hard-core bosons. Hole pairs, 
which appear with doping, have enough empty states to move, and they are responsible for 
low-energy kinetics. Above T c a material such as Y Ba2Cu^OQ +x contains a non-degenerate 
gas of hole bipolarons in singlet and in triplet states. Triplets are separated from singlets 
by a spin gap J and have a lower mass due to a lower binding energy, Fig. 12. The main 
part of the electron-electron correlation energy (Hubbard U and the inter-site Coulomb 
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FIG. 12: Bipolaron picture of high temperature superconductors. A corresponds to a singlet 
intersite bipolaron. B is a triplet intersite bipolaron, which naturally includes the addition of an 
extra excitation band. The crosses are copper sites and the circles are oxygen sites, w is a half 
bandwidth of the polaron band, t is a half bandwidth of the bipolaronic band, A/2 is the bipolaron 
binding energy per polaron and J is the exchange energy per bipolaron. 



repulsion) and the electron-phonon interaction are taken into account in the binding energy 
of bipolarons A, and in their band-width renormalization as described above. When the hole 
density is small, n& <C 1 (as in cuprates), their bipolaronic operators are almost bosonic. 
The hard-core interaction does not play any role in this dilute limit, so only the long- 
distance Coulomb repulsion is relevant. This repulsion is significantly reduced due to a large 
static dielectric constant in oxides, e ^> 1. Hence, carriers are almost free charged bosons 
coexisting with thermally excited nondegenerate single fermions, so that the conventional 
Boltzmann kinetics (see below) and the Bogoliubov transformation 58[ for a charged Bose 
gas are perfectly applied in the normal and in the superconducting state, respectively. 

The population of singlet, n s , triplet n t , and polaron, n p bands is determined by the 
chemical potential /i = Tiny, where y is found using the thermal equilibrium of singlet and 



triplet bipolarons and polarons, 

2n s + 2n t + n p = x. (57) 

Applying the effective-mass approximation for quasi-two-dimensional energy spectra of all 
particles we obtain for < y < 1 

-m* s * In (1 - y) - 3m* t * In (l - ye~ J/T ) + m* In (l + <//2 e -A/(2T)) = ™ (5g) 

in the normal state, and y = 1 in the superconducting state. Here x is the total number of 
holes per unit area. If the polaron energy spectrum is (quasi)one-dimensional, an additional 
T -1 / 2 appears in front of In in the third term on the left hand side of Eq.(59). 



D. Role of disorder and the phase diagram of cuprates 

We should also take into account localization of holes by the random potential, because 
doping inevitably introduces some disorder. The Coulomb repulsion restricts the number 
of charged bosons in each localized state, so that the distribution function will show a 
mobility edge E c j^. The number of bosons in a single potential well is determined by 
the competition between their long-range Coulomb repulsion c.a. 4e 2 /(eo£) and the binding 
energy E c — e. If the localization length diverges with the critical exponent v < 1 (£ ~ 
(E c — e)~ v ), we can apply a 'single well-single particle' approximation assuming that there 
is only one boson in each potential well. Within this approximation localized charged bosons 
effectively obey the Fermi-Dirac statistics, so that their density is given by 

n M T ) = / — i , 7 , (59) 



y- 1 exp(E/T) + 1 

where Nl(E) is the density of localized states. Near the mobility edge it remains constant 
Nl{E) ~ til/T, where T is of the order of the binding energy in a single random potential 
well, and ul is the number of localized states per unit area. The number of empty localized 
states turns out to be linear as a function of temperature in a wide temperature range T < T 
from Eq.(60). Then the conservation of the total number of carriers yields for the chemical 
potential: 

<% = -mT In (1 - y) - 3mf In (l - ye~ J / T ) + (60) 

m* In (1 + y V2 e -A/(2T)) _ 27mr ^ + ^ 
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FIG. 13: Phase diagram of superconducting cuprates in the bipolaron theory (courtesy of J. Hofer). 

If the number of localized states is about the same as the number of pairs, ~ x/2, 
a solution of this equation does not depend on temperature in a wide temperature range 
above T c . With y to be a constant (y ~ 0.6 in a wide range of parameters in Eq.(61)), the 
number of singlet bipolarons in the Bloch states is linear in temperature, 



The numbers of triplet pairs and single polarons are exponentially small at low temperatures, 



The model suggests a phase diagram of the cuprates as shown in Fig. 13. This phase 
diagram is based on the assumption that to account for the high values of T c in cuprates one 
has to consider electron-phonon interactions larger than those used in the BSC-Eliashberg 
theory of superconductivity. Regardless of the adiabatic ratio, the Migdal-Eliashberg theory 
of superconductivity and the Fermi-liquid theory break at A ~ 1. The many-electron system 
collapses into the small (bi)polaron regime at A > 1 with well separated vibrational and 
charge- carrier degrees of freedom. Even though it seems that these carriers should have 
a mass too large to be mobile, the inclusion of the on-site Coulomb repulsion and a poor 
screening of the long-range electron-phonon interaction leads to mobile intersite bipolarons 




(61) 



T < J, A/2. 



as discussed above. In the normal state the Bose gas of the bipolarons is non-degenerate 
and below T c their phase coherence sets in and hence superfluidity of the double-charged 2e 
bosons occurs. There are also thermally excited single polarons and triplets in the model, 
Fig. 12, which are responsible for the crossover regime at T* and normal state charge and spin 
(pseudo)gaps in cuprates. These pseudogaps were predicted as half of the binding energy A 



and the singlet-triplet exchange energy J of preformed bipolarons, respectively 
temperature, where the polaron density compares with the bipolaron one. 



23. T* is a 



E. Low Fermi energy: individual pairing in cuprates 

Clear evidence for bipolarons comes from a parameter-free estimate of the renormal- 
ized Fermi-energy ep, which yields a very small value in cuprates, where the band structure 
is quasi-two-dimensional with a few degenerate hole pockets. Applying the parabolic ap- 
proximation for the band dispersion we obtain the renormalized Fermi energy as 

£f = — ~, (62) 

m\ 

where d is the interplane distance, and nj,m* are the density of holes and their effective 
mass in each of the hole subbands i renormalized by the electron-phonon interaction. One 
can express the renormalized band-structure parameters through the in-plane magnetic-field 
penetration depth at T « 0, measured experimentally: 



As a result, we obtain a parameter-free expression for the "true" (i.e. renormalized) Fermi 
energy as 

£f = -: 2\2 ' ( 64 ) 
where g is the degeneracy of the spectrum, which may depend on doping in cuprates. One 
expects 4 hole pockets inside BZ due to the Mott-Hubbard gap in underdoped cuprates. If 
the hole band minima are shifted with doping to BZ boundaries, all their wave vectors would 
belong to the stars with two or more prongs. The groups of wave vectors for these stars 
have only ID representations. It means that the spectrum will be degenerate with respect 
to the number of prongs which the star has, i.e g > 2. The only exception is the minimum 
at k = (71", 7i") with one prong and g — 1. Hence, in cuprates the degeneracy is 1 < g < 4. 



Because Eq.(65) does not contain any other band-structure parameters, the estimate of ep 
using this equation does not depend very much on the parabolic approximation for the band 
dispersion. 

Generally, the ratios n/m* in Eq.(63) and in Eq.(64) are not necessarily the same. The 
'superfluid' density in Eq.(64) might be different from the total density of delocalized carriers 
in Eq.(63). However, in a translation invariant system they must be the same This 
is also true even in the extreme case of a pure two-dimensional superfluid, where quantum 
fluctuations are important. One can, however, obtain a reduced value of the zero temperature 
superfluid density in the dirty limit, I <C £(0), where £(0) is the zero-temperature coherence 
length. The latter was measured directly in cuprates as the size of the vortex core. It 
is about 10 A or even less. On the contrary, the mean free path was found surprisingly 
large at low temperatures, / ~ 100-1000 A. Hence, it is rather probable that all cuprate 
superconductors are in the clean limit, I ^> £(0), so the parameter-free expression for ep, 
Eq.(65), is perfectly applicable. 

A parameter-free estimate of the Fermi energy obtained by using Eq.(65) is presented 
in Table 1. The renormalized Fermi energy in many cuprates is less than 100 meV, if the 
degeneracy g > 2 is taken into account. That should be compared with the characteristic 
phonon frequency, which is estimated as the plasma frequency of oxygen ions, 



4irZ 2 e 2 N , . 

W = V^^' (65) 

One obtains LO=&AmeV with Z = 2, iV = 6/V ce ii, M = 16a.it. for YBCI2CU3OQ. Here V ce u 
is the volume of the (chemical) unit cell. The estimate agrees with the measured phonon 
spectra. As established experimentally in cuprates (see the introduction) high-frequency 
optical phonons are strongly coupled with holes. A low Fermi energy is a serious problem 
for the Migdal-Eliashberg approach. Since the Fermi energy is small and phonon frequencies 
are high, the Coulomb pseudopotential /x* is of the order of the bare Coulomb repulsion, 
fi* ~ /i c ~ 1 since the Tolmachev-Morel-Anderson logarithm is ineffective. Hence, to get a 
pairing, one has to have a strong coupling, A > fi c . However, one cannot increase A without 
accounting for the polaron collapse of the band. Even in the region of the applicability of 
the BCS-Eliashberg theory (i.e. at A < 0.5), the non-crossing diagrams cannot be treated as 
vertex corrections like in Ref. 62j ]. since they are comparable to the standard terms, when 
ojje F > 1. 



TABLE I: The Fermi energy (multiplied by the degeneracy) of cuprates 
Compound T c (K) \ H , a b (A) d(A) ge^ (meV) 
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In many cases (Table 1) the renormalized Fermi energy is so small that pairing is certainly 
individual, i.e. the bipolaron size is smaller than the inter-carrier distance. Indeed, this is 
the case, if 

e F < ttA. (66) 
If the bipolaron binding energy is twice of the pseudogap experimentally measured in the 



normal state of many cuprates 24j, A > lOOmeV, Eq.(67) is well satisfied in underdoped 
and even in a few optimally and overdoped cuprates. One should notice that the coherence 
length in a charged Bose gas has nothing to do with the size of the boson. It depends on 
the interparticle distance and the mean-free path 3], and might be as large as in the BCS 
superconductor. Hence, it would be incorrect to apply the ratio of the coherence length 
to the inter-carrier distance as a criterium of the BCS-Bose liquid crossover. The correct 
criterium is given by Eq.(67). 



IV. NORMAL STATE PROPERTIES OF CUPRATES IN FCM 



The low-energy FCM electronic structure of cuprates is shown in Fig. 12 |57[. Polaronic 
p-holes are bound in lattice inter-site singlets (A) or in singlets and triplets (B) (if spins 
are included in Eq.(2)) at any temperature. Above T c a charged bipolaronic Bose liquid is 
non-degenerate and below T c phase coherence (ODLRO) of the preformed bosons sets in. 
The state above T c is perfectly " normal" in the sense that the off-diagonal order parameter 
(i.e. the Bogoliubov-Gor'kov anomalous average jF(r, r') = (^(r)^(r')) is zero above the 
resistive transition temperature T c . Here ip^(r) annihilates electrons with spin j, j at point 
r. (Bi)polarons and thermally excited phonons are well decoupled in the strong-coupling 
regime of the electron-phonon interaction, so the conventional Boltzmann kinetics for mobile 
polaronic and bipolaronic carries is applied. 

A. Normal state in-plane resistivity, the Hall effect, magnetic susceptibility and 
the Lorenz number 

A nonlinear temperature dependence of the in-plane resistivity below T*, a temperature- 
dependent paramagnetic susceptibility, and a peculiar maximum in the Hall ratio well above 
T c have remained long-standing problems of cuprate physics. The bipolaron model pro- 



vides their quantitative description 



59 



63, 



64j ]. Here we use a 'minimum' bipolaron model 



Fig.l2A, which includes the singlet bipolaron band and the spin 1/2 polaron band sepa- 
rated by T*, and the r— approximation in weak electric E and magnetic fields, B _L E [64]. 
Bipolaron and single-polaron non-equilibrium distributions are found as 

/(k) = f (E) + r||v ■ {F + 0n x F} , (67) 



where v =dE/dk, F = V(/i - 2e</>), f (E) = [y- 1 exp(E/T) - and the Hall angle 
O = @{, = 2eBTb/rrib for bipolarons with the energy E = k 2 /(2rrib), and F = V(/i/2 — e<p), 
f (E) = {y-V 2 exp[(E + T*)/T\ + E = k 2 /(2m p ), and 6 = 6 P = eBr p /m p for ther- 

mally excited polarons. Here and m p are the bipolaron and polaron mass, respectively, 
y = exp(/i/T), fi is the chemical potential, and n = B/.B is a unit vector in the direction of 
the magnetic field. Eq.(68) is used to calculate the electrical resistivity and the Hall ratio 
as 

P Ae 2 r b n b (l + An p /n b y 1 ' 

1 + 2A 2 n p /n b 

H ' 2en b {l + An p /n b r 1 ' 

where A = T v rrib / (irbfrip) . The atomic densities of quasi two-dimensional carriers are found 

as 

n b = ^\Hl-y)\, (70) 
n =T£. In [1 + exp (_T7T)] . (71) 

7T 

and the chemical potential is determined by doping x using 2n& + n p = x — til, where hl 
is the number of carriers localized by disorder (here we take the lattice constant a = 1). 
Polarons are not degenerate. Their number remains small compared with twice the number 
of bipolarons, n p /{2rib) < 0.2, in the relevant temperature range T < T*, so that 

ywl-exp(-To/T), (72) 

where T = ir(x — n^/rrib ~ T c is about the superconducting critical temperature of the 
(quasi)two-dimensional Bose gas. Because of this reason, the experimental T c was taken as 
T in our fits. Using Eqs. (70,71,72) we obtain 

r (T\ — R i + 2Ay/ 2 (T/r e )exp(-rvr) 

Rh[T) - Rm {l + A(T/T c )y^e W (-T*/T)r (?3) 
where i?i/o — [ e ( x — n L)}^ 1 - If we assume that the number of localized carriers depends 
only weakly on temperature in underdoped cuprates since their average ionization energy 
is sufficiently large, then Rho is temperature independent at T < T*. As proposed in 



Ref. 



the scattering rate at relatively high temperatures is due to inelastic collisions of 



itinerant carriers with those localized by disorder, so it is proportional to T . We also 




have to take into account the residual scattering of polarons by optical phonons, so that 
t -1 oc (T/Ti) 2 + exp (— uj/T), if the temperature is low compared with the characteristic 
phonon energy u. The relaxation times of each type of carriers scales with their charge e* 
and mass as r Ptb oc m p 3 / 2 (e*)~ 2 , so we estimate A = (m b /m p ) 5 ^ 2 ps 6 if we take m b ps 2m p . 
As a result the in-plane resistivity is given by 

o(T) = o (T/T^ + eKpj-u/T) 

P{ ' P0 [l + A(T/T c )y 1 / 2 exp(-T*/T)Y { } 
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where po ~ mf,/[2e 2 (x — til)] and T\ are temperature independent. Finally, one can easi 
obtain the uniform magnetic susceptibility due to nondegenerate spin 1/2 polarons as [65] 

X (T)=V /2 exp(-T7T) 



Xo, 



(75) 



where B = (fi 2 B m p /7r), and xo is the magnetic susceptibility of the parent Mott insulator. 
Our model numerically fits the Hall ratio, Rh{T), the in-plane resistivity, p(T), and the 
magnetic susceptibility x{T) of Y Ba2CuzOT_$ within the physically relevant range of all 
parameters (see Fig. 14 and the Table). The ratio of polaron and bipolaron mobilities A = 7 
used in all fits is close to the above estimate, and y n ~ 1.5 x 10^ 4 emu/mole is very close to 
the susceptibility of a slightly doped insulator |66| |. The maximum of Rh(T) is due to the 
contribution of thermally excited polarons into transport, and the temperature dependence 
of the in-plane resistivity below T* is due to this contribution and the combination of the 
carrier-carrier and carrier-phonon scattering. The characteristic phonon frequency from the 
resistivity fit (see the Table) decreases with doping and the pseudogap T* shows the doping 



241- 



behavior as observed in other independent experiments 

Notwithstanding our explanation of the Hall ratio, the in-plane resistivity and the bulk 
magnetic susceptibility might be not so convincing as a direct measurement of the double 
charge 2e on carriers in the normal state. In 1993, we discussed the thermal conductivity of 
preformed bosons ((37 1 . The contribution from carriers to the thermal transport provided by 
the Wiedemann- Franz law depends strongly on the elementary charge as ~ (e*) -2 and should 



be significantly suppressed if e* = 2e. The Lorenz number, L, has been directly measured 
in YBa 2 Cu 3 OQ,g5 by Zhang et al. [68] using the thermal Hall conductivity. Remarkably the 
measured value of L just above T c was found just the same as predicted by the bipolaron 
model [gtI, L « 0.15L e , where L e is the conventional Fermi-liquid Lorenz number. The 
breakdown of the Wiedemann-Franz law has been also explained in the framework of the 

n 

bipolaron model [69j . 

B. Normal-state Nernst effect 

In disagreement with the weak-coupling BCS and the strong-coupling bipolaron theories 
a significant fraction of research in the field of high-temperature superconductivity suggests 
that the superconducting transition is only a phase ordering while the superconducting 
order parameter !F(r, r') remains nonzero above the resistive T c . One of the key experiments 
supporting this viewpoint is the l arg e Nernst signal observed in the normal (i.e. resistive) 
state of cuprates (see Ref. Q, E3 Q and references therein). Some authors Q 
claim that numerous resistive determinations of the upper critical field, H c2 (T) in cuprates 
have been misleading since the Nernst signal [3| and the diamagnetic magnetization [7^ 
imply that H c2 (T) remains large at T c and above. They propose a "vortex scenario", where 
the long-range phase coherence is destroyed by mobile vortices, but the amplitude of the 
off-diagonal order parameter remains finite and the Cooper pairing with a large binding 
energy exists well above T c supporting the so-called "preformed Cooper-pair" or "phase 
fluctuation" model [74]. The model is based on the assumption that the superfluid density is 
small compared with the normal carrier density in cuprates. These interpretations seriously 
undermine many theoretical and experimental works on superconducting cuprates, which 
consider the state above T c as perfectly normal with no off-diagonal order, either long or 
short. 

We believe that the vortex (or phase fluctuation) scenario contradicts straightforward 
resistive and other measurements, and it is theoretically inconsistent. This scenario is im- 
possible to reconcile with the extremely sharp resistive transitions at T c in high-quality 
underdoped, optimally doped and overdoped cuprates. For example, the in-plane and out- 
of-plane resistivity of Bi — 2212, where the anomalous Nernst signal has been measured 



is perfectly "normal" above T c , Fig. 15, showing only a few percent positive or negative 
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FIG. 15: In-plane (A) and out-of-plane (B) resistivity of 3 single crystals of Bi2Sr2CaCu20j. 
showing no signature of phase fluctuations above the resistive transition. 



□ 



magnetoresistance [j^l . explained with bipolarons [zf| • Both in-plane [zil S, 13, S l&j] 
and out-of-plane |82j, |S3> IS4 resistive transitions of high-quality samples are sharp and re- 
main sharp in the magnetic field providing a reliable determination of the genuine H c2 (T). 
The vortex entropy |7l| estimated from the Nernst signal is an order of magnitude smaller 
than the difference between the entropy of the superconducting state and the extrapolated 
entropy of the normal state obtained from the specific heat. The preformed Cooper-pair 
model [z4| is incompatible with a great number of thermodynamic, magnetic, and kinetic 
measurements, which show that only holes (density x), doped into a parent insulator are 
carriers both in the normal and the superconducting states of cuprates. The assumption 
[7^ that the superfluid density is small compared with the normal-state carrier density is 
also inconsistent with the theorem jfjl], which proves that the number of supercarriers at 
T = OK should be the same as the number of normal-state carriers in any clean superfluid. 



Recently we described t 
normal state phenomenon 



le^ unusual Nernst signal in cuprates in a different manner as the 
We have also extended this description to cuprates with very 



low doping level accounting for their Nernst signal, the thermopower and the insulating-like 



in-plane low temperature resistance 



861 ] as observed 



71 



72] 



Thermomagnetic effects appear in conductors subjected to a longitudinal temperature 
gradient V X T in x direction and a perpendicular magnetic field in z direction. The transverse 
Nernst-Ettingshausen effect 0] (here the Nernst effect) is the appearance of a transverse 
electric field E y in the third direction. When bipolarons are formed in the strong-coupling 



regime, the chemical potential is negative, Eq.(73). It is found in the impurity band just 
below the mobility edge at T > T c . Carriers, localized below the mobility edge contribute 
to the longitudinal transport together with the itinerant carriers in extended states above 
the mobility edge. Importantly the contribution of localized carriers of any statistics to the 
transverse transport is normally small [88] since a microscopic Hall voltage will only develop 
at junctions in the intersections of the percolation paths, and it is expected that these are 
few for the case of hopping conduction among disorder-localized states js^. Even if this 
contribution is not negligible, it adds to the contribution of itinerant carriers to produce 
a large Nernst signal, e y (T, B) = —E y fV x T, while it reduces the thermopower S and the 
Hall angle G. This unusual "symmetry breaking" is completely at variance with ordinary 
metals where the familiar " Sondheimer" cancelation j^] makes e y much smaller than 5* tan 
because of the electron-hole symmetry near the Fermi level. Such behavior originates in the 
"sign" (or "p — n") anomaly of the Hall conductivity of localized carriers. The sign of their 
Hall effect is often opposite to that of the thermopower as observed in many amorphous 
semiconductors 88] and described theoretically 

The Nernst signal is expressed in terms of the kinetic coefficients cr^- and as 

_ ^xx^yx Gyx^-xx (WRX 

e y~ a 2 +a 2 ' w 

v xx 1 xy 

where the current density is given by ji = cr^Ej + aijVjT. When the chemical potential p is 
at the mobility edge, localized carriers contribute to the transport, so and Oy in Eq.(77) 
can be expressed as a^ 1 + a\^ and a^ + cx ij, respectively. Since the Hall mobility of carriers 
localized below p, a l yx , has the sign opposite to that of carries in the extended states above 
p, (Ty^, the sign of the off-diagonal Peltier conductivity a l yx should be the same as the sign 
of a^f. Then neglecting the magneto-orbital effects in the resistivity (since 6 •€ 1 [ill) we 
obtain 

StanG « piatt - |«Ll)(©^ - I©'|) (77) 



® xx ® xy 



and 



e y ~P«* + Kl)-5tane, (78) 



where Q ext = cr^f /o~ xx , Q l = a l yx /a xx , and p = l/cr xx is the resistivity. 

Clearly the model, Eqs. (78,79) can account for a low value of StanO compared with a 

nn 

large value of e y in some underdoped cuprates [TiJ, 122] due to the sign anomaly. Even in 
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FIG. 16: Normal state in-plane resistivity of underdoped Lai.94 Sro.o6Cu04 (triangles [7l|]) as 
revealed in the field B = 12 Tesla and compared with the bipolaron theory, Eq.(82) (solid line). 



the case when localized bosons contribute little to the conductivity their contribution to 
the thermopower S = pjaj^ — \oi l xx \)) could almost cancel the opposite sign contribution 
of itinerant carriers [85(. Indeed the longitudinal conductivity of itinerant two-dimensional 
bosons, a ext oc J Q dEEdf(E)/dE diverges logarithmical when p in the Bose-Einstein distri- 
bution function f(E) = [exp((E — p)/T) — 1] _1 goes to zero and the relaxation time r is a 
constant. At the same time a e xx oc J Q dEEiE — p)df (E) / dE remains finite, and it could have 
the magnitude comparable with a l xx . Statistics of bipolarons gradually changes from Bose 
to Fermi statistics with lowering ener gy a cross the mobility edge because of the Coulomb 
repulsion of bosons in localized states 92]. Hence one can use the same expansion near the 
mobility edge as in ordinary amorphous semiconductors to obtain the familiar textbook re- 
sult S = SoT with a constant So at low temperatures |23|. The model becomes particularly 
simple, if we neglect the localized carrier contribution to p, and a xy , and take into account 
that a e X y oc B/p 2 and Q ext oc B/p in the Boltzmann theory. Then Eqs. (78,79) yield 

Stan6ocT/p (79) 

and 

e y (T,B)x(l-T/T 1 )/p. (80) 

A C c 0rdl „ g to our e arlier SUgg e Stl o n Q the insulating low-tomperature dependence of 
p(T) in underdoped cuprates originates from the elastic scattering of nondegenerate itinerant 
carriers off charged impurities. We assume here that the carrier density is temperature 



independent at low temperatures in agreement with the temperature-independent Hall effect 
95^ . The relaxation time of nondegenerate carriers depends on temperature as r oc T -1 / 2 
for scattering off short-range deep potential wells, and as T 1//2 for very shallow wells [94]. 
Combining both scattering rates we obtain 

p = p [(T/T 2 ) 1 / 2 + (T 2 /T) 1 / 2 ]. (81) 

Eq.(82) with p = 0.236 mO-cm and T 2 = 44. 6K fits extremely well the experimental 
insulating-like normal state resistivity of underdoped Lai.94 Sr .06CuO4 in the whole low- 
temperature range from 2K up to 50K, Fig. 16, as revealed in the field B = 12 Tesla j^l 3|- 
Another high quality fit can be obtained combining the Brooks-Herring formula for the 



3D scattering off screened charged impurities, as proposed in Ref. 



for almost undoped 



LSCO, or the Coulomb scattering in 2D (r oc T) and a temperature independent scattering 



rate off neutral impurities with the carrier exchange [97[ similar to the scattering of slow 
electrons by hydrogen atoms in three dimensions. Hence the scale T 2 , which determines the 
crossover toward an insulating behavior, depends on the relative strength of two scattering 
mechanisms. Importantly our expressions (80,81) for StanO and e y do not depend on the 
particular scattering mechanism. Taking into account the excellent fit of Eq.(82) to the 
experiment, they can be parameterized as 



Stan 9 = e ^ T ^,^ , (82) 
1 + T/T 2 ' v ; 

and 

e y {T, B) = e t 2 + T ' ^ ' 

where T\ and eo are temperature independent. 

In spite of all simplifications, the model describes remarkably well both S"tanB and 
e y measured in Lai.94 Sr .06CuO 4 with a single fitting parameter, T\ = 50K using the 
experimental p(T). The constant eo = 2.95 pV/K scales the magnitudes of StanO and 
e y . The magnetic field B — 12 Tesla destroys the superconducting state of the low-doped 
La x 94 Sr 06 CuO 4 down to 2K, Fig. 16, so any residual superconducting order above 2K is 
clearly ruled out, while the Nernst signal, Fig. 17, is remarkably large. The coexistence of the 
large Nernst signal and a nonmetallic resistivity is in sharp disagreement with the vortex 
scenario, but in agreement with our model. Taking into account the field dependence of 
the conductivity of localized carriers, the phonon-drug effect, and their contribution to the 



1 3/2 




transverse magnetotransport can well describe the magnetic field dependence of the Nernst 
signal and improve the fit in Fig. 17 at the expense of the increasing number of fitting 
parameters. 



C. Normal state diamagnetism 

A number of experiments (see, for example, 000 [l0 Eml, and references 
therein), including torque magnetometries, showed enhanced diamagnetism above T c , which 
has been explained as t he fl uctuation diamagnetism in quasi-2D superconducting cuprates 



(see, for example Ref. 



10l(). The data taken at relatively low magnetic fields (typically 



below 5 Tesla) revealed a crossing point in the magnetization M (T, B) of most anisotropic 
cuprates (e.g. Bi - 2212), or in M(T,B)/B^ 2 of less anisotropic YBCO 0. The depen- 
dence of magnetization (or M/B 1 ^ 2 ) on the magnetic field has been shown to vanish at some 
characteristic temperature below T c . However the data taken in high magnetic fields (up to 
30 Tesla) have shown that the crossing point, anticipated for low- dimensional superconduc- 
tors and associated with superconducting fluctuations, does not explicitly exist in magnetic 
fields above 5 Tesla 

Most surprisingly the torque magnetometery H uncovered a diamagnetic signal 



somewhat above T c which increases in magnitude with applied magnetic field. It has been 
linked with the Nernst signal and mobile vortexes in the normal state of cuprates [73^ . 
However, apart from the inconsistences mentioned above, the vortex scenario of the normal- 
state diamagnetism is internally inconsistent. Accepting the vortex scenario and fitting the 
magnetization data in Bi — 2212 with the conventional logarithmic field dependence Q|, one 
obtains surprisingly high upper critical fields > 120 Tesla and a very large Ginzburg- 
Landau parameter, k — A/£ > 450 even at temperatures close to T c . The in-plane low- 
temperature mag netic field penetration depth is A = 200 nm in optimally doped Bi — 2212 
(see, for example |l03^ ). Hence the zero temperature coherence length £ turns out to be about 
the lattice constant, £ = 0.45nm, or even smaller. Such a small coherence length rules out 
the "preformed Cooper pairs" [74j, since the pairs are virtually not overlapped at any size of 
the Fermi surface in Bi — 2212 . Moreover the magnetic field dependence of M(T, B) at and 
above T c is entirely inconsistent with what one expects from a vortex liquid. While — M(B) 
decreases logarithmical at temperatures well below T c , the experimental curves 0, 
clearly show that — M(B) increases with the field at and above T c , just opposite to what 
one could expect in the vortex liquid. This significant departure from the London liquid 



behavior c 
transition 



early indicates that the vortex liquid does not appear above the resistive phase 



SSI- 



Some time ago we explained the anomalous diam agnetism in cuprates as the Landau 
normal-state diamagnetis m of preformed bosons loJ |. The same model predicted the un- 
usual u pper critical field |l05 | observed in many superconducting cuprates 0, 0, s, 
Isil . Is^ Il0f| (see below). Here we extend the model to high magnetic fields taking into 
account the magnetic pair-breaking of singlet bipolarons and the anisotropy of the energy 
spectrum. When the strong magnetic field is applied perpendicular to the copper-oxygen 
plains the quasi-2D bipolaron energy spectrum is quantized as 

E a = u(n + 1/2) + 2t c [l - cos{K z d)}, (84) 

where u = 2e5/m{,, n = 0, 1,2, and t c , K z , d are the hopping integral, the momentum 
and the lattice period perpendicular to the planes. Quantum numbers a also include the 
momentum along one of the in-plane directions. Expanding the Bose-Einstein distribution 
function in powers of exp[([i — E a )/T] with the negative /1 one can readily obtain (after 



summation over n) the boson density 

eB 

and the magnetization 



iid 1 — exp(— uik/T ) 



(85) 



M(T, B) 



eT 
ltd 



-'"-/"< + — ^fofock/T) — 



exp[(/2 - 2t c )fc/T] 



fc=i 



x 



(86) 



exp(— uok/T) 

1 a; exp(— uk/T) \ 
k ~ T[l -exp{-uk/T)]) ' 

where /if, = e/m^, /2 = u —u/2 and /oO^) is the modified Bessel function. At low temperatures 
T — > Schafroth's result 107| is recovered, M(0, B) = — n&/i&. The magnetization of charged 
bosons is field-independent at low temperatures. At high temperatures, T ^> T c the chemical 
potential has a large magnitude , so we can keep only terms with k — 1 in Eqs. (86,87) to 
obtain 



M(T, B) 



-ribUb + 



Trih 



1 



u;exp(— lo/T) 



(87) 



£ V' T[l-esxp(-u/T)] i 
The experimental conditions are such that T ^> uj when T is of the order of T c or higher, so 

that 

M(T,B) = -n bf i b ^, (88) 

which is the Landau orbital diamagnetism of nondegenerate carriers. The bipolaron in-plane 
mass in cuprates is about m& ~ 10m e [3]. Using this mass yields M(0, B) w 2000 A/m with 
the bipolaron density nf, = 10 21 cm -3 . Then the magnitude and the field/temperature 
dependence of M (T, B) near and above T c are about the same as experimentally observed 



in Refs 



73 



The pseudogap temperature T* depends on the magnetic field predominantly because of 
the magnetic-field splitting of the single-polaron band in Fig. 12. As a result the bipolaron 
density depends on the field (as well as on temperature) near T c as 



n b (T, B) = n b (T c , 0) 1 + (T e - T)/f - (B/B ) 



(89) 



where To and Bo are constants depending on T*, (3 = 2 if the polaron spectrum is spin- 
degenerate, and (3 = 1 if the spin degeneracy is removed by the crystal field already in the 
absence of the external field. 

Theoretical temperature and field depende ncie s of M(T,B), Eq.(87) agree qualitatively 



with the experimental curves in Bi — 2212 [73, 



1001 ] . if the depletion of the bipolaron density, 



Eq.(90) is taken into account. The depletion of rib accounts for the absence of the crossing 
point in M(T, B) at high magnetic fields. Nevertheless a quantitative fit to experimental 
M(T, B) curves using T and B as the fitting parameters is premature. The experimen- 
tal diamagnetic magnetization has been extracted from the total magnetization assuming 
that the normal state paramagnetic contribution remains temperature- independent at all 



temperatures 



100]. This assumption is inconsistent with a great number of NMR and 
the Knight shift measurements, and even with the preformed Cooper-pair model itself. 
The Pauli spin-susceptibility has been found temperature-dependent in these experiments 
revealing normal-state pseudo gaps, contrary to the assumption. Hence the experimental 
diamagnetic M(T, B) has to be corrected by taking into account the temperature 

dependence of the spin paramagnetism at relatively low temperatures. 



D. Spin pseudogap, c-axis transport and charge pseudogap 

The pairing of holes into singlets well above T c should be seen as a drop of the nuclear 
magnetic relaxation rate 1/Ti with temperature lowering. Indeed it is a common feature 
of the normal state of many cuprates. The bipolaron model has described the temperature 
dependence of 1/T\ |23[. The conventional contact hyperfine coupling of nuclear spin on a 
site i with electron spins is given in the site representation by 

H i = A^c^ + H.c, (90) 

3 

where Ai is an operator acting on the nuclear spin, and j is its nearest neighbor sites. 
Performing projecting transformations to bipolarons as above we obtain the effective spin- 
flip interaction of triplet bipolarons with the nuclear spin as 

H * ~ E b ) M + H - c - ( 91 ) 

Here 1,1' = 0, ±1 are z-components of spin S — 1. The NMR width due to the spin-flip 
scattering of triplet bipolarons on nuclei is obtained using the Fermi-Dirac golden rule, 

1 B f 2t , df(E) , . 

where f(E) = [exp((E + J — fj)/T)— ] _1 is the triplet distribution function, and 2t is their 
bandwidth. For simplicity the triplet DOS is taken as a constant (= l/(2i)). As a result, 



we obtain 



BTsmh(t/T) 



(93) 



T x t 2 [cosh[(i + J)/T - lny] - cosh(t/T)] ' 
where B is a temperature independent hyperfine coupling constant. 

Eq.(94) describes all essential features of the nuclear spin relaxation rate in copper-based 
oxides: the absence of the Hebel-Slichter coherent peak below T c , the temperature-dependent 
Korringa ratio (1/TTi) above T c , and a large value of 1/ Tj d ue to the small bandwidth 
21 It nicely fits the experimental data in Y Ba 2 Cu<±0% |l08 | with reasonable values of 
the parameters, t = 250K and J = 150K 2^. A similar unusual behavior of NMR was 
found in underdoped Y Ba 2 Cu 3 0^ +x ^ and in many other cuprates. The Knight shift, which 
measures the spin susceptibility of carriers, also drops well above T c in many copper oxides, 
in agreement with the bipolaron model . Th e 'spin' gap has been observed above and below 
T c in Y Ba 2 Cu ?J 0§ +x with unpolarized 109| . and polarized 110 1 neutron scattering. 



The bipolaron model 



of cuprates 



83 



84 



ras also quantitatively explained c-axis transport and the anisotropy 



The crucial point is that single polarons dominate in c-axis 



transport at finite temperatures because they are much lighter than bipolarons in c-direction. 
Bipolarons can propagate across the planes due to a simultaneous two-particle tunnelling 
alone, which is much less probable than a single polaron tunnelling. Along the planes 
polarons and inter-site bipolarons propagate with comparable effective masses, as shown 
above. Hence in the mixture of nondegenerate quasi-two-dimensional bosons and thermally 
excited fermions, only fermions contribute to c -axis transport, if the temperature is not 
very low, which leads to thermally activated c -axis transport and to a fundamental relation 
between the anisotropy and the uniform magnetic susceptibility of cuprates |25| . 

The exponential temperature dependence of c-axis resistivity and "c" versus " ab" 
anisotropy was interpreted within t he fr amework of the bipol aron model in many cuprates, 



in particular in La 2 _ x Sr x CuO/ X (25, 



113 . Bi 2 Sr 2 CaCu 2 



8+5 



mj, YBa 2 Cu 3 6+x [84], and 



HgBa 2 CuO±+s [83]. Importantly, the uniform magnetic susceptibility above T c increases 
with doping. It proves once more that cuprates are doped insulators, where low energy 
charge and spin degrees of freedom are due to holes doped into a parent insulating ma- 
trix with no free carriers and no free spins. A rather low magnetic susceptibility of parent 
insulators in their paramagnetic phase is presumably due to a singlet pairing of copper spins. 



V. SUPERCONDUCTING STATE OF CUPRATES 



A. Parameter-free evaluation of T c : Bose-Einstein condensation versus the 
Kosterlitz-Thouless transition 



An ultimate goal of the theory of superconductivity is to provide an expression for T c as 
a function of some well-defined parameters characterizing the material. In the framework 
of the BCS theory T c is fairly approximated by the familiar McMillan's formula, which 
works well for simple metals and their alloys. But applying a theory of this kind to high-T c 
cuprates is problematic. Since bare electron bands are narrow, strong correlations result 
in the Mott insulating state of undoped parent compounds. As a result, fi* is ill-defined 
in doped cuprates, and polaronic effects are important as in many doped semiconductors. 
Hence, an estimate of T c in cuprates within the BCS theory appears to be an exercise in 
calculating the Coulomb pseudopotential /i* rather than T c itself. One cannot increase A 
either without accounting for a polaron collapse of the band as discussed above. This appears 
at A rs 1. 

On the other hand, the bipolaron theory provides a parameter-free expression for T< 



m, 



which fits the experimentally measured T c in many cuprates for any level of doping. T c is 
calculated using the density sum rule as the Bose-Einstein condensation (BEC) temperature 
of 2e charged bosons on a lattice. Just before the discovery |2fi | we predic ted T r as high as 
w 100K using an estimate of the bipolaron effective mass |l!4 |. Uemura |l lfs| established 
a correlation of T c with the in-plane magnetic field penetration depth measured by fisR 
technique in many cuprates as T c ~ 1/A^ 6 . The technique is based on the implantation of 
spin polarized muons. It monitors the time evolution of the muon spin polarization. He 
concluded that cuprates are neither BCS nor BEC superfluids but they are in a crossover 
region from one to the other, because the experimental T c was found about 3 or more times 
below the BEC temperature. 

Here we calculate T c of a bipolaronic superconductor taking properly into account the 
microscopic band structure of bipolarons in layered cuprates as derived in s ection 3. We 



115j involves not 



arrive at a parameter-free expression for T c , which in contrast to Ref. 
only the in-plane, A a j, but also the out-of-plane, A c , magnetic field penetration depth, and 
a normal state Hall ratio Rh just above the transition. It describes the experimental data 
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FIG. 18: Theoretical critical temperature compared with the experiment (the theory is ex- 
act for samples on the straight line) for LaSrCuO compounds (squares), for Zn substituted 
Y Ba2Cui- x Zn x 07 (circles), for Y Ba2CuzO'j_$ (triangles), and for HgBa2CuOi + s (diamonds). 
Experimental data for the London penetration depth are taken from T. Xiang et al, Int. J. 
Mod. Phys. B12, 1007 (1998) and B. Janossy et al, Physica C181, 51 (1991) in YBa 2 Cu 3 7 s 
and Y Ba2Cu\- x Zn x O-j; from V.G. Grebennik et al, Hyperfine Interactions 61, 1093 (1990) and 
C. Panagopoulos (private communication) in underdoped and overdoped La2- x Sr x CuO^, respec- 
tively, and from J. Hofer et al, Physica C, 297, 103 (1998) in HgBa2Cu04 + $. The Hall coefficient 
above T c is taken from A. Carrington et al, Phys. Rev. B48, 13051 (1993) and J. R. Cooper 
(private communication) (Y ' Ba 2 Cu^O-j^s an d Y Ba2Cu\- x Zn x O'j) and from H.Y. Hwang et al, 
Phys. Rev. Lett. 72, 2636 (1994) (La 2 - x Sr x CuO A ). 




for a few dozen different samples, Fig. 18, clearly indicating that many cuprates are in the 
BEC rather than in the crossover regime. The energy spectrum of bipolarons is at least 
two- fold degenerate in cuprates (section 3). One can apply the effective mass approximation 
at T ~ T c , Eq.(52), because T c should be less than the bipolaron bandwidth. Also three- 
dimensional corrections to the spectrum are important for the Bose-Einstein condensation. 
They are well described by the tight-binding approximation as 

E * = + i^f + 2tc[1 - cos(i ^ )] - (94) 

Substituting the spectrum, Eq.(95) into the density sum rule, 

[eME l K /T c ) - If 1 = n b (95) 

K,«=(x,2/) 
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FIG. 19: Correction coefficient to the 3D Bose-Einstein condensation temperature as a function of 
anisotropy. 



one readily obtains T c as (in ordinary units) 



k B T c = f 



x 



3.31£ 2 (n 6 /2) 2/3 



(96) 



k B T c J {m**m**m**) 1 /^ ' 

where the coefficient f(x) ps 1 is shown in Fig. 19 as a function of the anisotropy, t c /(k B T c ) : 
and m** = /l 2 /(2|t c |d 2 ). 

This expression is rather ambiguous because the effective mass tensor as well as the 
bipolaron density rib are not well known. Fortunately, we can express the band-structure 
parameters via in-plane, 

A a 6 



m* x *m* y * 



87m;,e 2 (m** + m**) 



1/2 



and out-of-plane penetration depths, 



A, 



m. 



1/2 



lQirnbe 2 

(we use c = 1). The bipolaron density is expressed through the in-plane Hall ratio (above 
the transition) as 



R 



H 



X 



4m* x *m** 



which leads to 



2enb {m** + m**) 2 ' 
T c = 1.64/ 



/,. \ / eR H x 1/3 



k B T c 



KM 



(97) 



(98) 



Here T c is measured in Kelvin, eRu in cm 3 and A in cm. The coefficient / is about unity 
in a very wide range of t c /(k B T c ) > 0.01, Fig.19. Hence, the bipolaron theory yields a 



parameter-free expression, which unambiguously tells us how near cuprates are to the BEC 
regime, 

/ p D \l/3 

T c « T C (3D) = 1.64 . (99) 

We compare two last expressions with the experimental T c of more than 30 different cuprates, 
for which both \ a b and A c are measured along with Rh(T c + 0) in Table 2 and Fig. 18. The Hall 
ratio has a strong temperature dependence above T c . Therefore, we use the experimental Hall ratio 
just above the transition. In a few cases (mercury compounds), where Rh(T c + 0) is unknown, 
we take the inverse chemical density of carriers (divided by e) as Rh- For almost all samples the 
theoretical T c fits experimental values within an experimental error bar for the penetration depth 
(about ±10%). There are a few Zn doped YBCO samples, Fig. 18, whose critical temperature is 
higher than the theoretical estimate. If we assume that the degeneracy of the bipolaron spectrum 
is removed by the random potential of Zn, then the theoretical T c would be almost the same as 
the experimental values for these samples as well. 

One can argue that due to a large anisotropy cuprates may belong to a 2D 'JCY' u niver s ality 
class with the Kosterlitz-Thouless (KT) critical temperature Tkt of preformed bosons |ll(J . IhtI | 



or the Cooper pairs 



74| . Should it be the case, one would hardly discriminate the Cooper pairs 



with respect to bipolarons. KT critical temperature is expressed through the in-plane penetration 



' , P 

depth alone as 7J] 



It appears significantly higher than the experimental values in many cases (see Table 2). There 
are also quite a few samples with about the same X a b and the same d, but with very different 
values of T c , which proves that the phase transition is not the KT transition. On the contrary, our 
parameter-free fit of the experimental critical temperature and the critical behavior (see below) 
favor 3D Bose-Einstein condensation of charged bosons as the mech anism of high T c rather than any 
low-dimensional phase-fluctuation scenario. The fluctuation theory jl isl ] further confirms the three- 
dimensional character of the phase transition in cuprates. However, it does not mean that all high- 
temperature superconductors are in the BEC regime with charged bosons as supercarriers. Some of 
them, in particular, electron-doped cuprates, MgBi and doped fullerenes might be in the BCS or 
intermediate regime, which makes the BCS-BEC cross over problem to be relevant. Starting from 
the pioneering works by Eagles 



cross 

3 



119] and Legget |l20l | this problem received particular attention 



TABLE II: Experimental data on T C (K), ab and c penetration depth(nm), Hall coeffi- 
cient (10 -3 (cm 3 /C)), and calculated values of T c respectively for La2- x Sr x CuO A (La), 
YBaCuO{x%Zn) (Zn), YBa 2 Cu 3 7 - x (Y) and HgBa 2 CuO A+x (Hg) compounds 

Compound T™ p X ab A c R H , T C (3D) T c T KT 
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TABLE III: Mass enhancement in cuprates 
Compound m c 



La(0.2) 


22.1 3558 


La(0.15) 


15.0 2698 


La(O.l) 


11.3 1909 


y(o.o) 


7.2 584 


F(0.12) 


8.3 600 


y(o.3) 


10.6 1994 



in the framew ork o f a negative Hubbard U model 

n 



path integral 



121 



122]. Both analytical (diagrammatic 



123], 



1241 ] ) and numerical 112511 studies have addressed the i ntermedi ate coupling regime 



125 



126, 



. However, in using 



beyond a variational approximation 121], including 2D systems 
the negative Hubbard U model, we have to realize that this model, which predicts a smooth BCS- 
BEC crossover, cannot be applied to the BCS-bipolaron crossover. The essential effect of the 
polaron band-narrowing (section 2) is missing in the negative (and positive) Hubbard U model. 
The polaron collapse of the bandwidth is mainly responsible for high T c . It strongly affects the 
BCS-BEC crossover significantly reducing the crossover region. 

It is interesting to estimate the effective mass tensor using the penetration depth and the Hall 
ratio. These estimates for in-plane and out-of-plane boson masses are presented in Table 3. They 
agree with the inter-site bipolaron mass (section 3). We notice, however, that the absolute value 
of the effective mass in terms of the free electron mass does not describe the actual band mass 
renormalization if the bare (band) mass is unknown. Nevertheless an assumption that the 
number of carriers is determined by the Luttinger theorem (i.e. n ~ 1) would lead to much heavier 
carriers with m* about 100m e . 



B. Isotope effect on T c and on supercarrier mass 

The advances in the fabrication of the isotope substituted samples made it possible to measure a 
sizable isotope effect , a = —d\aT c / d\a M in many high-T c oxides. This led to a general conclusion 
that phonons are relevant for high T c . Moreover the isotope effect in cuprates was found to be 
quite different from the BCS prediction, a = 0.5 (or less). Several compounds showed a > 0.5 



1281 ]. and a small negative value of a was found in Bi — 2223 

Essential features of the isotope effec t, in particular large values in low T c cuprates, an overall 



trend to lower value as T c increases 130], and a small or eve n neg ative a in some high T c cuprates 
were understood in the framework of the bipolaron theory [l3l] . With increasing ion mass the 
bipolaron mass increases and the Bose-Einstein condensation temperature T c decreases in the 
bipolaronic superconductor. On the contrary in polaronic superconductors an increase of the ion 
mass leads to a band narrowing and to an enhancement of the polaron density of states, and 
to an increase of T c . Hence the isotope exponent in T c can distinguish the BCS like polaronic 
superconductivity with a < , and the Bose-Einstein condensation of small bipolarons with a > 0. 
Moreover, underdoped cuprates, which are definitely in the BEC regime, could have a > 0.5, as 
observed. 



The isotope effect on T c is linked with the isotope effect on the carrier mass, a m * , as 



131] 



a = -dhxT c /d\o.M = a m *[l - Z/(X - /u c )], (101) 

where a m * = dlnm* /din Al and Z = m/m* <C 1. In ordinary metals, where the Migdal approx- 
imation is believed to be valid, the renormalized effective mass of electrons is independent of the 
ion mass M because the electron-phonon interaction constant A does not depend on M. However, 
when the e-ph interaction is sufficiently strong, the electrons form polarons dressed by lattice dis- 
tortions, with an effective mass m* = mexp^Ep/fkv). While E p in the above expression does not 
depend on the ion mass, the phonon frequency does. As a result, t here is a large isotope effect 
on the carrier mass in polaronic conductors, a. m * = (l/2)ln(m*/m) [laj . in contrast to the zero 
isotope effect in ordinary metals. 

Such an effect was observed in cuprates in the London penetration depth of isotope-substituted 
samples d2j. The carrier density is unchanged with the isotope substitution of O 16 by O 18 , so that 
the isotope effect on \ a b measures directly the isotope effect on the carrier mass. In particular, 
the carrier mass isotope exponent a m * was found as large as a m * = 0.8 in Lai.895>S'ro.i05Cu04. 
Then the polaron mass enhancement should be m**/m ~ 5 in this material. Using Eq.(57) 
we obtain the in-plane bipolaron mass as large as m** ~ 10m e with the bare hopping integral 
T(NNN) = 0.2 eV. The in-plane magnetic field penetration depth, calculated with this mass 
is X a b = [m**/87rne 2 ] 1 / 2 ~ 316nm, where n is the hole density. It agrees well with the experi- 
mental one, X a b ~ 320nm. Using the measured values of X a b = 320 nm, A c = 4160 nm, and of 



TABLE IV: Coherence volume Q in A 3 , the in-plane £ a 6 and out- of- plane £ c coherence lengths 
derived from a Ginzburg-Landau analysis of the specific heat 

Compound tt £^,(A 2 ) £ c , (A) 

YBa 2 Cu 3 O r 400 125 3.2 

YBa 2 Cu 3 7 -o. 2 5 3 09 119 2.6 

FSa2C-ii30 7 _o.o5 250 119 2.1 

YBa 2 Cu 3 O 7 - .i 143 119 1.2 

Cao.8Fo.2Sr 2 T/ .5^o.5C«2C>7 84 70 1.2 

Th. 8 Ba2Ca 2 .2Cu 3 O 10 40 < 0.9 

= 4 x 10~ 3 cm? /C (just above T c ) we obtain T c = 2>1K from Eq.(100) in astonishing agree- 
ment with the experimental value T c = 30 K in this compound. More recent high resolution angle 
resolved photoemission spectroscopy studies [l3] provided further compelling evidence for strong 
e-ph interaction in the cuprates. They revealed a fine phonon structure in the electron self-energy 
of underdoped La2_ x Sr x .Cu04 samples and a complicated isotope effect in the electron spectral 
function of Bi2212 that depended on the electron energy and momentum. 



C. Specific heat anomaly 



Bose liquids (or more precisely He 4 ) show the characteristic A -point singularity of their specific 
heat, but superfluid Fermi liquids like BCS superconductors exhibit a sharp second order phase 
tran s ition accompanied by a finite jump in the specific heat. It was established beyond doubt 



132 



133, 



134, 



135, 



136] that the anomaly in hig h cuprates differs qualitatively from the BSC 
prediction. As was stressed by Salamon et al. |l37l ] the heat capacity is logarithmic near the 
transition, and consequently, cannot be adequately treated by the mean-field BCS theory even 
including the gaussian fluctuations. In particular, estimates using the gaussian fluctua tion s yield 



an unusually small coherence volume, Table 4 



magnetic field de pendence of the anomaly 



bipolaron model 



106 



comparable with the unit cell volume 



1331] . The 



138] is also unusual, but it can be described by the 



Calculations of the specific heat of charged bosons in a magnetic 
field require an analytical DOS, N(e,B) of a particle, scattered by other particles and/or by a 
random potential of impurities. One can use DOS in the magnetic field with an impurity scattering 
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FIG. 20: Temperature dependence of the specific heat divided by temperature (arb. units) of 
the charged Bose-gas scattered by impurities for several fields {ujh = 2eB/m**). (b) shows two 
anomalies, the lowest one traces resistive transition, while the highest anomaly is the normal state 
feature. 



calculated in the non-crossing approximation 

C(T,B) _ d 



1051 ] . The specific heat coefficient 
de- N ^ B > 



(102) 



T TdT J ~~exp[(e-//)/T] - 1' 

calculated with this DOS and with ii determined from rib = f deN(e, B)f(e), is shown in Fig. 20. 
The broad maximum at T ~ T c is practically the same as in the ideal Bose gas without scattering. 
It barely shifts in the magnetic field. However, there is the second anomaly at lower temperatures, 
which is absent in the ideal gas. It shifts with the magnetic field, tracing precisely the resistive 
transition, as clearly seen from the difference between the specific heat in the field and zero-field 
curve, Fig. 20b. The specific heat, Fig. 20, is in st rikin g resemblance with the Geneva group's 
experiments on DyBa2Cu^0i and on Y Ba2Cu^O-j where both anomalies were observed. 

Within the bipolaron model, when the magnetic field is applied, it hardly changes the temperature 
dependence of the chemical potential near the zero field T c because the energy spectrum of thermally 
excited bosons is practically unchanged. This is because their characteristic energy (of the order 
of T c ) remains huge compared with the magnetic energy of the order of 2eB/m** . In contrast, the 
energy spectrum of low energy bosons is strongly perturbed even by a weak magnetic field. As a 
result the chemical potential 'touches' the band edge at lower temperatures, while having almost 
the same 'kink'-like temperature dependence around T c as in zero field. While the lower anomaly 



corresponds to the true long-range order due to the Bose-Einstein condensation, the higher one 
is just a 'memory' about the zero-field transition. This microscopic consideration shows that a 
genuine phase transition into a superconducting state is related to a resistive transition, and to 
the lower specific heat anomaly, while the broad higher anomaly is a normal state feature of the 
bosonic system in the external magnetic field. Different from the BCS superconductor these two 
anomalies are well separated in the bosonic superconductor at any field but zero. 



D. Universal upper critical field 



The upper critical field, H C 2(T) = <I>o/27r£(T) 2 , is very different in the BCS superconductor 
and in the charged Bose-gas (CBG). While H C 2(T) is linear in temperature near T c in the Landau 
theo ry of second-order phase transitions, it has a positive curvature H C 2(T) ~ (T c — T) 3 / 2 in CBG 



105]. Also at zero temperature H C 2(0) is normally below the Pauli pair-breaking limit given by 
Hp ~ 1.84T C (in Te sla) in the BCS the ory, but the limit c an b e exceeded by many times in C BG. 



In cuprates 



77 



78 



79, 



80 



81 



85, 



1061 ] . spin-ladders 



1401 ] and organic superconductors |l41 | 



high magnetic field studies revealed a non-BCS upward curvature of resistive H r <?j T). When mea- 



78, 



79. (1, Q[l41|, the 



surements were performed on low-T c unconventional superconductors 
Pauli limit was exceeded by several times. A non-linear temper ature dependence in the vicinity 

106]. Importantly, a thermody- 



80 



81 



85 



1421 ] due to contrasting 



of T c was unambiguously observed in a few samples 
namically determined H C 2 turned out much higher than the resistive H C 2 
magnetic field dependencies of the specific heat anomaly and of resistive transition. 

I believe that many unconventional superconductors are in the 'bosonic' limit of preformed real- 
space bipolarons, s o the ir resistive H C 2 is actually a critical field of the Bose-Einstein condensation 



of charged bosons 1051 ] . Calculations above carried out for the heat capacity of CBG lead to the 
conclusion that the resistive H C 2 and the thermodynamically determined H C 2 are very different 
in bosonic superconductors. While the magnetic field destroys the condensate of ideal bosons, it 
hardly shifts the specific heat anomaly as observed. 

A comprehen sive scaling of resistive H C 2 measurements in unconventional superconductors is 
shown in Fig. 21 106J] in the framework of the microscopic model of charged bosons scattered by 
impurities. An expression for H C 2(T) accounting for a temperature dependence of the number of 




FIG. 21: Resistive upper critical field (determined at 50% of the trans ition') of cuprates, spin- 
ladders and organic superconductors scaled according to Eq.(105) 10 J |. The parameter b is 1 
(solid line), 0.02 (dashed-dotted line), 0.0012 (dotted line), and (dashed line). The inset shows a 
universal scaling of the same data near T c on the logarithmic scale. 



delocalized bosons, nj,(T), can be written as 

H c2 (T) = H 



106] 



n b (T) 



,1/2 



3/2 



(103) 



tn b (T c ) 

where T c is the zero-field critical temperature, and t = T/T c . Here the scaling constant Hq 
depends on the mean-free path I , Hq = ^>q/2ttS,q, with the characteristic (coherence) length 
£o — GA^C^c)) 1 ■ I n the vicinity of T c one obtains the parameter-free H C 2(T) oc (1 — i) 3 / 2 using 
this equation, but the low-temperature behavior depends on a particular scattering mechanism, 
and a detailed structure of the density of localized states. As suggested by the normal state Hall 
measurements in cuprates n&(T) can be parameterized as n&(T) = nb(0) + constant x T (see also 
Eq.(62)), so that H C 2{T) is described by a single-parameter expression as 



H c2 (T) = H 



6(1 -t) 



+ l-tV2 



3/2 



(104) 



The parameter b is proportional to the number of delocalized bosons at zero temperature. We 
expect that this expression is applied in the whole temperature region except ultra-low tempera- 



143]. Exceeding the Pauli pair-breaking limit 



tures, where the non-crossing approximation fails 
readily follows from the fact, that the singlet-pair binding energy is related to the normal-state 
pseudogap temperature T* , rather than to T c . T* is higher than T c in bosonic superconductors 
and in cuprates. The universal scaling of H C 2 near T c is confirmed by resistive measurements of 
the upper critical field of many cuprates, spin-ladders, and organic superconductors, as shown in 
Fig. 21. All measurements reveal a universal (1 — 1) 3//2 behavior in a wide temperature region (inset) 
near T c . The low-temperature behavior of H C 2{T) / Hq is not universal, but well described using 
Eq.(105) with the single fitting parameter, b. The parameter is close to 1 in high quality cuprates 
with a very narrow resistive transition. It naturally becomes rather small in overdoped cuprates 
where randomness is more essential. 



E. Symmetry and space modulations of the order parameter 



Independent observations of normal state pseudogaps in a number of magnetic and kinetic 
measurements, and the unusual critical behavior, discussed above, tell us that many cuprates may 
not be BCS superconductors. Indeed their superconducting state is as anomalous as the normal 
one. In particular, there is strong evidence for ad-like order parameter (changing sign when 
the CuO? plane is rotated by tt/2) in cuprates A number of phase-sensitive experiments 



145] provide una mbiguous evidence in this direction; furthermore, the low temperature magnetic 



penetration depth 



146, 



1471 ] was found to be linear in a few cuprates as expected for a d-wave BCS 



superconductor. However, SIN and SIS tunnelli ng s tudies, the c-axis Josephson tunnelling 
and some high-precision magnetic measurements |l43 | show a more usual s-like symmetry or even 
reveal an upturn in the temperature dependence of the penetration de pth b elow some characteristic 



temperature 



|l5fj. 



Also both angle-resolved photoemission (ARPES) 



microscopy (STM) 



151] and scanning tunnelling 



1521 ] have shown that the maximum energy gap and 2A/T C ratio is several times 



larger than expected in the weak-coupling BCS theory, or in its intermediate-coupling Eliashberg 
generalization. Strong deviations from the Fermi/BCS-liquid behavior are suggestive of a new 
electronic state in cuprates, which is a charged Bose liquid of bipolarons 

Actually there are more complicated deviations from the conventional Fermi/BCS-liquid be- 
havior than the normal state pseudogaps. Recent studies of the gap function revealed two dis- 
tinctly different gaps with different magnetic field and temperature dependence and the 



156] a nd w ithout 



157 



15a ] applied 



153 ] explaining two different 



checkerboard spatial modulations of the tunnelling DOS, with 
magnetic fields. We have proposed a simple phenomenological model 
gaps in cuprates. The main assumption, supported by a parameter-free estimate of the Fermi 
energy (section 3.5), is that the attractive potential is large compared with the renormalized Fermi 
energy, so that the ground state is the Bose-Einstein condensate of tightly bound real-space pairs. 
Here I present an e xpla nation of the symmetry of the order parameter and real-space modulations 



of tunnelling DOS 



1601 ] in the framework of the bipolaron theory. 



The anomalous Bogoliubov-Gor'kov average J-(r,r') depends on the relative coordinate of two 
electrons (holes) and on the center-of-mass coordinate, p = r — r',R = (r + r')/2. Its Fourier 
transform, /(k, K), depends on the relative momentum k and on the center-of-mass momentum 
K. In the BCS theory K =0, and the Fourier transform of the order parameter is proportional 
to the gap in the quasi-particle excitation spectrum, /(k, K) ~ Ak. Hence the symmetry of the 
order parameter and the symmetry of the gap are the same in the weak-coupling regime. Under 
the rotation of the coordinate system, A^ changes its sign, if the Cooper pairing appears in the 
d-channel. 

Real-space pairs might also have an unconventional symmetry due to a specific symmetry of 
the pairing potential as in the case of the Cooper pairs, but in any case the ground state and 
DOS are homogeneous, if pairs are condensed with K = 0. On the other hand, the symmetry 
of the order parameter could be different from the 'internal' symmetry of the pair wave function, 
and from the symmetry of a single-particle excitation gap in the strong-coupling regime [!| • If the 
pair band dispersion has its minima at finite K in the center-of-mass BZ, the Bose condensate is 
inhomogeneous. In particular, the center-of-mass bipolaron energy bands could have the ir m inima 
at the Brillouin zone boundaries at K = (tt, 0) and three other equivalent momenta 160] (the 



lattice constant is taken as a = 1). These four states are degenerate, so that the condensate wave 
function ?/>( m ) i n the real (Wannier) space, m = (m x ,m y ), is their superposition, 

V>(m)= £ b K e~ lKm , (105) 

K=(±7r,0),(0,±7r) 

where &k = ±\/wc/2 are c-numbers, and n c (T) is the atomic density of the Bose-condensate. The 
superposition, Eq.(106), respects the time-reversal and parity symmetries, if 

Tp(m) = ■sfn c [cos(7rm x ) ± cos(-7rm y )] . (106) 

The order parameter, Eq.(107), has d-wave symmetry changing sign in the real space, when the 



lattice is rotated by tt/2. This symmetry is entirely due to the pair-band energy dispersion with 
four minima at K / 0, rather than due a specific pairing potential. It reveals itself as a checker- 
board modulation of the hole density with two-dimensional patterns, oriented along the diagonals. 
From this i nsigh t one can expect a fundamental connection between st ripe s detected by different 



techniques 



160] 



1611 ] and the symmetry of the order parameter in cuprates 
Importantly, even if preformed singlet bipolarons are condensed at T point of their center-of- 
mass BZ (i.e. with K = 0), the superconducting order-parameter could be d — w ave due to the 



1621 ] the orientation 



"orientation" degeneracy of inter-site pairs, Fig. 10. As proposed by Andreev 
degeneracy of inter-site pairs provides d-wave symmetry also irrespective of the symmetry of the 
pairing potential, but without any space modulations of the hole density. It could be relevant for 
overdoped cuprates (see below). 

Now let us take into account that in the superconducting state (T < T c ) single-particle exci- 
tations interact with the pair condensate via the same attractive potential, which forms the pairs 



159]. The Hamiltonian describing the interaction of excitations with the pair Bose-condensate in 



the Wannier representation is 

H = -J2 [*( m " n ) + Mm,n]4mC S n + ^[A(m)c{ m q m + H.c], (107) 

s,m,n m 

where s =T, I is the spin, and A(m) oc ip(m) is a coherent gap function. Applying equations of 
motion for the Heisenberg operators cj m (t) and c sm (t), and the Bogoliubov transformation 



.(*) = ^K(m)a,e- fc - t + u^m)/^*"*], (108) 



q m (i) = £>„(m)A,e-^* - <(m)at e ^*], (i 09 ) 
one obtains BdG equations describing the single-particle excitation spectrum, 

eu{ra) = - ^[t(m - n) + ^ m , n ]u(n) + A(m>(m), (110) 

n 

-ev(m) = - ^J[i(m - n) + n5 m ,n]v(n) + A(m)«(m), (111) 
n 

where excitation quantum numbers v are omitted for transparency. Different from the conven- 
tional BdG equations in the weak-coupling limit, there is virtually no feedback of single particle 
excitations on the off-diagonal potential, A(m), in the strong-coupling regime. The number of 



these excitations is low at temperatures below T* = A p , so that the coherent potential A(m) is an 



159]. 



external (rather than a self-consistent) field, solely determined by the pair Bose condensate 
While the analytical solution is not possible for any arbitrary off-diagonal interaction A(m), one 
can readily solve the infinite system of discrete equations (111,112) for a periodic A(m) with a 
period commensurate with the lattice constant. For example, 

A(m) = A c [e inm * - e i7rw »], (112) 

corresponds to the pair condensate at K = (±7r,0) and (0, ±7r), Eq.(107), with a temperature 



dependent (coherent) A c oc WnjT). In this case the quasi-momentum k is the proper quantum 
number, v = k, and the excitation wave-function is a superposition of plane waves, 

u u (m) = u k e ik ' m + Uke^-e)" 1 , (113) 
^(m) = V]{ e i(k ~ s ^ m + v ic e i(k - s y> m . (114) 

Here g x = (tt, 0), g y = (0,7r), and g = (7r,7r) are reciprocal doubled lattice vectors. Substituting 
these equations into the BdG equations (111,112) one obtains four coupled algebraic equations, 

ek^k = ^k - A c (-u k - v k ), (115) 

£k«k = £k- g ^k + A c (v k - v k ), (116) 

-ek^k = 6c-g^k + A c (n k - tt k ), (117) 

-ek^k = Ck-g^k - A c (u k - tt k ), (118) 

where £ k = — X^ n ^( n ) e * k n — A 4 - The determinant of the system (116-119) yields the following 
equation for the energy spectrum e: 



(e ~ £k)(e - &- g )(e + £k- g J(e + £k- g 



y ' 



= A^2e + £ k _ gx +£ k _ gj/ )(2 e -£ k -£ k _ g ). (119) 

Two positive roots for e describe the single-particle excitation spectrum. Their calculation is rather 

cumbersome, but not in the extreme strong-coupling limit, where the pair binding energy 2A p is 

large compared with A c and with the polaron bandwidth. The chemical potential in this limit is 

pinned below a single-particle band edge, so \i is negative, and its magnitude is large compared 

with A c . Then the right hand side in Eq.(120) is a perturbation, and the spectrum is 

A 2 

eik ~ £k - — , (120) 
A 2 

e 2 k ~ £k- g - — • (121) 



If a metallic tip is placed at the point m above the surface of a sample, the STM current 
I(V, m) creates an electron (or hole) at this point. Applying the Fermi-Dirac golden rule and the 
Bogoliubov transformation, Eqs.(109,110), and assuming that the temperature is much lower than 
Ap/kB one readily obtains the tunnelling conductance 

a(V, m) = dI{ ^ v m) oc K(m)\ 2 HeV - e„), (122) 
which is a local excitation DOS. The solution Eq.(114) leads to a spatially modulated conductance, 

a(V, m) = a reg (V) + a mod {V) cos(-Km x + irm y ). (123) 
The smooth (regular) contribution is 

<Jre g (V)=(T (Urk + Urk)S(eV - € rk ), (124) 

k,r=l,2 

and the amplitude of the modulated contribution is 

cr mod (V)=2a ^ u rk u rk 5(eV - e rk ), (125) 
k,r=l,2 

where do is a constant. Conductance modulations reveal a checkerboard pattern, as the Bose 
condensate itself, 

— — = A cos(irm x + Trm y ) , (126) 

&reg 

where 

A = 2^ [u lk u lk S(eV - eik) + M 2 k« 2 k<5(eF - e 2 k)] / 

k 

+ «ik)<*(eV - ei k ) + (n| k + <4)<*(eV - e 2k )] 

k 

is the amplitude of modulations depending on the voltage V and temperature. An analytical result 

can be obtained in the strong-coupling limit with the excitation spectrum given by Eqs. (121,122) 

for the voltage near the threshold, eV « A p . In this case only states near bottoms of each excitation 

band contribute to the integrals in Eq.(125), so that 

£k - eik A 2 

ui k = uik « -uik < «ik, (127) 

eik - Ck-g 



and 



£k- g -e 2k _ _ A c _ 
u 2k = 5 — — U2k ~ _ U2k ^ U2k . (128) 

e 2 k - £k fJ-w 



Substituting these expressions into A, Eq.(127), yields in the lowest order of A c , 

2A? 



11.111 



(129) 



156, 



157, 



158 



163], where spatial checkerboard 



The result, Eq.(127), is reminiscent of STM data 
modulations of a were observed in a few cuprates. Both commensurate and incommensurate 
modulations were found depending on sample composition. In our model the period is determined 
by the center-of mass wave vectors K of the Bose-condensed preformed pairs. While the general 
case has to be solved numerically, the perturbation result, Eq.(127) is qualitatively applied for 
any K at least close to T c , where the coherent gap is small, if one replaces cos^m^ + Trm y ) by 
cos(K x m x + KyUiy). Different from any other scenario, proposed so far, the hole density, which is 
about twice of the condensate density at low temperatures, is spatially modulated with the period 
determined by the inverse wave vectors corresponding to the center-of-mass pair ba nd-m inima. This 



'kinetic' interpretation of charge modulations in cuprates was originally proposed 160] before STM 
results became available. It could account for those DOS modulations in superconducting samples, 
which disappear above T c because the coherent gap A C (T) vanishes, so that ^4 = above T c in 
Eq.(127). Indeed some inelastic neutron scattering experiments show that in com mensurate inelastic 



peaks are observed only in the superconducting state of high-T c cuprates 



1641 ] . The vanishing at 



T c of incommensurate peaks is inconsistent with any other stripe picture, where a characteristic 



distance needs to 
(see, for example 



3e observed in the normal state as well. On the other hand some STM studies 



1651 ]) report incommensurate and com men surate DOS modulations somewhat 



1661 ]. I believe that those modulations 



above T c , in particular, in heavily underdoped cuprates 
are due to a single-particle band structure and impurity states near the top of the valence band in 
doped charge-transfer insulators, rather than a signature of any cooperative phenomenon. 

In this way the strong-coupling Frohlich-Coulomb model, Eq.(2), links charge heterogeneity, 
pairing, and pseudo-gaps as manifestations of the strong electron-phonon attractive interaction in 
narrow bands of doped Mott-Hubbard insulators. 



VI. OVERDOPED CUPRATES: BOSON-FERMION MIXTURES 

It has been mentioned in section 1.1 that the chemical potential could enter the oxygen band 
in overdoped samples 3] as a result of the overlap of bipolaron and polaron bands, so the Fermi- 
level crossing could be seen in ARPES, and pseudogaps gradually vanish, Fig.22. If the number of 
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FIG. 22: The Pauli exclusion principle prevents a decay of mobile bipolarons into single polarons 
in overdoped cuprates because all states below the Fermi level \i are occupied. 



states below the Fermi level in the polaron band is less than the number of holes, bipolarons remain 
stable mobile quasi-particles because the Pauli exclusion prin ciple prevents their decay into single 



1671 ] of a two-body problem on the 



polarons, Fig. 22. A comprehensive analysis by Kornilovitch 
square lattice with the nearest-neighbor attraction strongly supports such a model of overdoped 
cuprates. It has been found that the stability of pairs increases with their momentum. The pairs 
are formed easier along the (tt, 0) direction than along the (ir, tt) direction. This might lead to 
the appearance of "hot pairing spots" on the K x and K y axes, while other regions of BZ remain 
unpaired. 

Such a system represents a mixture of mobile 2e-charged bosons and e-charged degenerate 



168], somewhat similar to neutral He 3 -He 4 



fermions, first considered by us before the discovery 
mixtures. While the normal state kinetic properties are dominated by a lighter polaronic Fermi 
liquid, the superconducting state is the Bose-Einstein condensate (BEC) of heavier bipolarons in 
the mixture. As we discuss in the remaining part, it is valid even in the extreme case of intrinsically 
immobile bipolarons, if they are hybridized with mobile polarons, irrespective to the strength of 
the hybridization interaction. 



A. Mobile fermions hybridized with immobile bosons: Boson-Fermion model 



A s an alte r native to the picture o f HTS with mobile bipolarons, Figs. 12,22, some authors 



169 



170 



171 



172 



gnnn 



17.1 11741 . 117,4 Il7fl| proposed a so-called Boson-Fermion model (BFM), where 



intrinsically immobile pairs are coupled with the Fermi sea of itinerant single fermions by the 
hybridization interaction. Soon after Anderson and Street and Mott 53j introduced localized pairs 



in amorphous semiconductors, such two component model of negative U centers coupled with 
the Fermi sea of itine rant ferm ions was employed to study superconductivity in disordered metal- 



semiconductor alloys 



177, 



1781 ]. When the attractive potential U is large, the model is reduced to 



localized hard-core bosons spontaneously decaying into itinerant electro ns and vi ce versa, different 



162, 



168 



from a non-converting mixture of mobile charged bosons and fermions 

This boson- fermion model (BFM) was applied more generally to describe pairing electron pro- 



cesses with localization-delocalization 
high-t emperatur e superconductors 



Refs. 



174, 



3 



17S |. The mode l attracted more attention in c onnection with 



170 



171 



172 



173 



174 



175 



176, 



180 



1811 ] . In particular, 



1751 ] claimed that 2D BFM with immobile hard-core bosons is capable to reproduce 



some physical properties and the phase diagram of cuprates. BFM has been also ad opted for a 



182]. 



description of superfluidity of atomic fermions scattered into bound (molecular) states 

Most studies of BFM below its transition into a low-temperature condensed phase applied a 
mean- field approximation (MFA), replacing zero-momen tum boson operators by c-numbers and 



169, 



170, 



172 



iratorsDv c-nu 

i7al I17I [17I 



176, 



182]. 



neglecting the boson self-energy in the density sum rule 
When the bare boson energy is well above the chem ical potential, the BCS ground state was found 

MFA led to a conclusion that BFM exhibits 



with bosons being only virtually excited 



169, 



170 



features compatible with BCS characteristics 

n 

17q . The transition was found more mean-field-like than the usua 



172|, and describes a crossover from the BCS-like 

Bose 

1 va lue of the fluctuation parameter Gi 



173] 



to a local pair behavior 

condensation, i.e. characterized by a relatively sma 

At the same time our previous study of BFM 
the boso n se lf-energy on the normal state boson spectral function and the transition temperature 
T r . Ref. 



1801 ] beyond MFA revealed a crucial effect of 



180] proved that the Cooper pairing of fermions via virtual bosonic states is impossible 
in any-dimensional BFM. It occurs only simultaneously with the Bose-Einstein condensation of 
real bosons. The origin of this simultaneous condensation lies in a softening of the boson mode 
at T = T c caused by its hybridization with fermions. The energy of zero-momentum bosons is 



renormalized down to zero at T = T r . no matter how weak the boson- fermion coupling and how 
large the bare boson energy are 180J. One can also expect that the boson self-energy should 
qualitatively modify the phase transition and the condensed phase of BFM below T c . 



B. Absence of BCS-BEC crossover in BFM 

Let us examine the phase transition and the condensed state of BFM beyond the ordinary 
mean-field approximation in two (2D) and three (3D)dimensions. It appears that T c = K in the 
two-dimensional model, even in the absence of any Coulomb repulsion, and the phase transition 
is never a BCS-like second-order phase transition even in 3D BFM because of the complete boson 
softening. 

The 2D BFM is defined by the Hamiltonian, 

H = E £k4 iff c ki(T + £o5>^ q + (130) 

giV _1 / 2 ^ (Vk^C-k+qATCk+q^.l + ff.Cj , 
q,k 

where £k = — 2t(cosk x + cosk y ) — /x is the 2D energy spectrum of fermions, Eq = Ab — 2/j, is 
the bare boson energy with respect to their chemical potential 2/i, g is the magnitude of the 
anisotropic hybridization interaction, 0k = 0_k is the anisotropy factor, and N is the number of 
cells. Here and further I take h = c = ks = 1 and the lattice constant a = 1. Ref. 174| argued 



that 'superconductivity is induced in this model from the anisotropic charge-exchange interaction 
(g</>k) between the conduction-band fermions and the immobile hard-core bosons', and 'the on-site 
Coulomb repulsion competes with this pairing' reducing the critical temperature T c less than by 



25%. Also it has been argued 175J, that the calculated upper critical field of the model fits well 
the experimental results in cuprates. 

Here I show that T c = K in the two-dimensional model, Eq.(131), even in the absence of any 
Coulomb repulsion, and the mean-field approximation is meaningless for any-dimensional BFM 
because of the complete boson softening. Replacing boson operators by c-numbers for q = 
in Eq.(131) one obtains a linearized BCS-like equation for the fermion order-parameter (the gap 
function) Ak, 

Ak = ^V0k Ak/tanh( ^ /2r) , (131) 



with the coupling constant g 2 = g 2 (l — 2n B ), renormalized by the hard-core effects. Using a two- 
particle fermion vertex part in the Cooper channel one can prove that this equation is perfectly 
correct even beyond the conventional non-crossing approximation 180(. The problem with MFA 



does not stem from this BSC-like equation, but from an incorrect definition of the bare boson 
energy with respect to the chemical pot entia l, Eq(T). This energy is determined by the atomic 



density of bosons (rib) a s (Eq.(9) in Ref. 



m) 



tanh — - = 1 — 2nh. 
2T 



(132) 



While Eq.(132) is perfectly correct, Eq.(133) is incorrect because the boson self-energy X;,(q, 0) 



due to the same hybridization interaction is missing. At first sight 



1741 ] the self-energy is small in 



comparison to the ki netic energy of fermions, if g is small. However S;,(0, 0) diverges logarithmical 
at zero temperature 18o| . no matter how week the interaction is. Therefore it should be kept in 
the density sum-rule, Eq.(133). Introducing the boson Green's function 

1 - 2n b 



£>(q,fi) 

one must replace incorrect Eq.(133) by 



iQ — Eq — Sft(q, O) 



(133) 



-^e <nT Z>(q,fi)=n 6 , 



(134) 



q,n 



where r = +0, and = 2-kTu (n = 0, ±1, ±2...). 

The divergent (cooperon) contribution to £&(q, f2) is given by 



180], 



S 6 (q,J2) 



2N 



(135) 



tanh[£ k _ q/2 /(2T)] + tanh[e k+q/2 /(2T)] 



^k-q/2 + £k+q/2 ~ *H 



so that one obtains 



S 6 (q,0) =S b (0,0) + 



2M* 



+ 0(q 4 ) 



(136) 



for small q and any anisotropy factor compatible with t he po int-group symmetry of the cuprates. 



Here M* is the boson mass, calculated analytically in Ref. 



and parabolic fermion band dispersion (see also Ref. 
nontrivial solution for A k at T = T c , if 



180] for the isotropic exchange interaction 



18ll[). The BCS-like equation (132) has a 



£b = -E 6 (0,0). 



(137) 



Substituting Eqs. (137,138) into the sum-rule, Eq.(135), one obtains a logarithmical divergent in- 
tegral with respect to q, and 

m const 

T c = T -— = 0. 138 
Jo d l/l 

The devastating result, Eq.(139) is a direct consequence of the well-known the orem , whi ch states 
that BEC is impossible in 2D. One may erroneously believe that MFA results are still 

applied in three-dimensions, where BEC is possible. However, increasing dimensionality does not 
make MFA a meaningful approximation for the boson-fermion model. This approximation leads 



to a conclusion that a BCS-like superconducting state occurs below the critical temper ature T r ^ 



169 



/iexp {—Eq/zc) via fermion pairs being virtually excited into unoccupied bosonic states 
Here z c = g 2 N(Q) and N(0) is the density of states (DOS) in the fermionic band near the Fermi 
level /i. However, the Cooper pairing of fermions is impossible via virtual unoccupied bosonic states 
also in 3D BFM. Indeed, Eqs. (132, 138) do not depend on the dimensionality, so that the analytical 
continuation of Eq.(132) to real frequencies u> yields the partial boson DOS as p(uj) = (l — 2nb)S(uj) 
at T = T c and q = in any-dimensional BFM. The Cooper pairing may occur only sim ultan eously 
with the Bose-Einstein condensation of real bosons in the exact theory of 3D BFM 1801 ] . The 
origin of the simultaneous condensation of the fermionic and bosonic fields in 3D BFM lies in the 
softening of the boson mode at T = T c caused by its hybridization with fermions. 

Taking into account the boson damping and dispersion shows that the boson spectrum signifi- 
cantly changes for all momenta. Continuing the self-ene rgy, Eq . ( 136) to real frequencies yields the 



damping (i.e. the imaginary part of the self-energy) as 



180] 



7(q ' w) = 4^ ln 



cosh(g£ + u/(4T c )) 



(139) 



_cosh(-g£ + w/(4T c ))_ 

where £ = t>i?/(4T c ) is a coherence length, and vf is the Fermi velocity. The damping is significant 
when q£ « 1. In this region 7(q, to) = ujttz c /(8T c ) is comparable or even larger than the boson 
energy u. Henc e bosons look like overdamped diffusive modes, rath er th an quasiparticles in the 
long-wave limit contrary to the erroneous 

of coherent free-particle-like motion of the bosons' in this limit. Only outside the long-wave region, 
the damping becomes small. Indeed, using Eq.(138) one obtains 7(q, u) = ujirz c /(2qvF) « 
so that b oson s at q >> l/£ are well defined quasiparticles with a logarithmic dispersion, ui(q) = 
z c ln(g£) 18fl| . Hence the boson energy disperses over the whole energy interval from zero up to 

Eq. 



The main mathematical problem with MFA in 3D also stems from the density sum rule, Eq.(135) 
which determines the chemical potential of the system and consequently the bare boson energy 
Eq(T) as a function of temperature. In the framework of MFA one takes th e bar e boson energy 
in Eq.(132) as a temperature independent parameter, Eq = g 2 N(0)ln(fj,/T c ) or determines 



it from t 



170 



re cqnsery ation of the total number of particles neglecting the boson self-energy, Eq.(133) 



174 



equation 



176 



182j |. Then Eq.(132) looks like the conventional linearized Ginzburg-Landau-Gor'kov 
1831 ] with a negative coefficient a oc T — T c at T < T c in the linear term. Then one 
concludes that the pha se transition is almost the conventional BCS-like transition, at least at 
E » T c 



169, 



170, 



]. These findings are mathematically and physically erroneous. Indeed, the 
term of the sum in Eq.(135) with £l n = is given by the integral 

^ 1 (140) 



' 2tt3£ + S fe (q,0) 
The integral converges, if and only if Eq > — Sft(0, 0). In fact, 

£ + £ 6 (0,0)=0 (141) 

is strictly zero in the Bose-condensed state, because m = —[Eq + £&(0,0)] corresponds to the 
boson chemical potential relative to the lower edge of the boson energy spectrum. More generally, 
Hb = corresponds to the appearance of the Bogoliubov-Goldstone mode due to a broken symmetry 
below T c . This exact result makes the BSC equation (132) simply an identity with a(T) = at 
any temperature below T c . On the other hand, MFA violates the density sum-rule, predicting 
the wrong negative a{T) below T c . Since a{T) = 0, one may expect that the conventional upper 
critical field, H C 2(T) is zero in BFM. To dete rmin e H C 2(T) and explore the condensed phase of 3D 



BFM, one can apply the Gor'kov formalism 



183], as described below. 



C. Normal and anomalous Green's functions of 3D BFM: pairing of bosons 



Let us now explore a simpli fied v ersion of 3D BFM in an external magnetic field B = V x A 



neglecting the hard-core effects 



mi, 



H = dr^4(r)h(r)Mr)+g[Hr)i>\(r)i>\(r) + H.c.] 

J s 

+ £o^(r)<Hr), (142) 

where ^ s (r) and (j)(r) are fermionic and bosonic fields, s =|,| is the spin, h(r) = — [V + 
ieA(r)] 2 /(2m) — fx is the fermion kinetic energy operator. Here the volume of the system is taken 



asV = 1. 

The Matsubara field operators, Q = exp(iir)Q(r) exp(— Ht), Q = ex.p(Hr)Q' (r) exp(—Hr) 
(Q = ip s , <j)) evolve with the imaginary time —1/T < r < 1/T as 



^■•I (r - r) = ^(r.rl + ^r^M, (143) 



dr 

dr 
d<f)(r,T) 



h*(r)^(r,r) - 9 ^t,t)^(t,t), (144) 
£ <Kr,r) + #|(r,T)V-t(r,T). (145) 



The theory of the condensed state can be formulated with the normal and anomalous fermion GFs 



lSa, S(r,r',r) = -(T T ^(r,r)^ a (r' J 0)) J ^+(r,r',r) = (T T ^ i (r,r)^ r (r' J 0)>, respectively, where 



the operation T T performs the time ordering. Fermionic and bosonic fields condense simultaneously 



180]. Following Bogoliubov [58( the bosonic condensate is described by separating a large matrix 
element <^o( r ) hi 4>( r 7 T ) as a number, while the remaining part 4>(r, r) describes a supracondensate 
field, 0(r, r) = <po(r) + 0(r, r). Then using Eq.(146) one obtains 

#o(r) = A(r) = -fW, r, 0+), (146) 

where ^, T ) = (^(r^.O), _The ec.ua.ions f„ r GFs are obtained o y n sl n g E q , (U4- 
146) and the diagrammatic technique [18a ] in the framework of the non-crossing approximation 



1861 ] . as shown in Fig. 23 and Fig. 24. 

An important novel feature of BFM is a pairing of supracondensate bosons, caused by their 
hybridization with the fermionic condensate, as follows from the last diagram in Fig. 24. Hence, 
one has to introduce an anomalous supracondensate boson GF, B + (r, r',r) = (T T (fi(r,T)(J)(r' , 0)} 
along with the normal boson GF, T>(r, r',r) = — (T T <p(r, r)0(r / ,O)). 

The diagrams, Fig. 23 and Fig. 24, are transformed into analytical equations for the time Fourier- 
components of the fermion GFs with the Matsubara frequencies lo = nT(2n + l) (n = 0, ±1, ±2, ...) 



as 



[ico - h(r)]g u (r,T') = 5(r - r') - A(r)^+(r,r') 

LU — UJ 

cj(x, r ) 

g 2 Tj2 J dx^+(r,x)^ w (r,x)^+(x,r / ), (147) 



[-iu - h*(r)]F+(r,r') = A*(r)&,(r,r') 



< <^ As *- + < ; < f ^ 



FIG. 23: Diagrams for the normal and anomalous fermion GFs. Zig-zag arrows represent the 
single-particle Bose condensate 4>o, dotted lines are the boson GFs, solid lines are the fermion GFs. 
Vertex (dot) corresponds to the hybridization interaction. 



g 2 T^2 J dx^(r,x)P^_ tJ (x,r)^+(x,r / ) 

ui' 

g 2 T^2 J ^_ a; ,(r,x)Bl a; _ w/ (r,x)g a; (x,r / ), 

U)' 

(iQ-E )V n (r,r') = 6(r-r') 
g 2 T^J dx^(r,x)g^(r,x)P n (x,r / ) 

9 2 TY,J ^(r,x)/n- w '(r I x)B+(x,r'), (148) 



and 



+ 



(-ifi-£*>)B+(r,r') 



S^E / ^ w ,(r,x)^ n+ „,(r,x)2?n(x,r') 
~ 9 2 TY,J dx^(x,r)^_ n (x,r)£+(x,r'). 
for the boson GFs with 6(r,r',r) = (T r 0(r,r)0(r / , 0)). 



> •> => ■■>■ +> |>.....-...„ 



o 



+ > 4 ><■ •• 



■< •> = ■■■<■ • # ■> 

5+ 



o 



+ < 4 ><■ =>■ 



FIG. 24: Diagrams for the supracondensate boson GFs. The Cooper-pairing of fermions leads to 
the Cooper-pair-like boson condensate, described by the boson anomalous GF, B + . 

D. Gor'kov expansion 

These equations can be formally solved in the homogeneous case without the external field, 
A = 0. Transforming into the momentum space yields GFs' time-space Fourier components as 

=~ p w > (149) 

**™ - F H^Sk ^ (150) 

and 

p( q>w ) = __ s ^_+^° , (151) 

|^-E | 2 + |r(q,^)|2' 

= - , (152) 

|^-^ | 2 + |r(q,^)| 2 ' 

where a) = w+zS/(k, u;), = fi+i£{,(q, fi), and £k = k 2 /(2m) — \i. The fermionic order parameter, 
renormalized with respect to the mean-field A due to the formation of the boson-pair condensate, 
is given by 

A(k,u,) = A + 5 2 T^ J |^jr+(k- q ,u/)tf(q,u, + u/), (153) 



and the boson-pair order parameter, generated by the hybridization with the fermion Cooper pairs, 



is 



r(q,ft) = 5 2 T^ I ^(k > w , )J r (q-k,n-a/). (154) 

Hence, there are three coupled condensates in the model described by the off-diagonal fields g4>o, 
A, and T, rather than two, as in MFA. At low temperatures all of them have about the same 
magnitude, as the fermion and boson self-energies, 



Z f Qt,u;) = -g 2 Tj2 [ 0<?(q-k,-u/)2?(q,u,-u/), 

Ul' 

S 6 (q,0) = - 5 2 T^ f ^0(k,a/)0(q-k,n-a; / ), 



(155) 



(156) 



respectively. 

On the other hand, when the temperature is close to T c (i.e. T c — T <C T c ), the boson pair 
condensate is weak compared with the single-particle boson and the Cooper pair condensates. In 
this temperature range T, Eq.(155) is of the second order in A, T oc A 2 , so that the anomalous 
boson GF can be neglected, since A is small. The fermion self-energy, Eq.(156) is a regular function 
of u and k, so it can be absorbed in the renormalized fermion band dispersion. Then the fermion 



normal and anomalous GFs, Eqs. 
can apply the Gor'kov expansion 



150 ,151) look like the familiar GFs of the BCS theory, and one 



1831 ] in powers of A(r) to describe the condensed phase of BFM 
in the magnetic field near the transition. Using Eq.(147) one obtains to the terms linear in A 

A*(r) = £r]T J dx^l(x,r)A*(x)gW(x,r). (157) 

The spatial variations of the vector potential are small near the transition. If A(r) varies slowly, 
the no rmal state GF, Q!^\r,r') differs from the zero-field normal state GF, Q^\r — r') only by a 



phase 



1831 ] (r, r') = exp[— zeA(r) • (r — r')]^ (r — r'). Expanding all quantities near the point 
x = r in Eq.(158) up to the second order in x — r inclusive, one obtains the linearized equation for 
the fermionic order parameter as 

7[V - 2ieA(r)] 2 A(r) = «A(r), (158) 



where 



£/Q -C/Q 1 



and 7 « 7C(3)t> 2 ,5 2 Ar(0)/(487r 2 T 2 £ 



The coefficient a{T) disappears in Eq.(159), since Eq = — £&(0, 0) at and below T c . It means that 
the phase transition is never a BCS-like second-order phase transition even at large Eq and small g. 
In fact, the transition is driven by the Bose-Einstein condensation of real bosons with q = 0, which 
occur due to the complete softening of their spectrum at T c . Remarkably, the conventional upper 
critical field, determined as the field, where a non-trivial solution of the linearized Gor'kov equation 



174 ] using MFA. 



(159) occurs, is zero in BFM, H C 2{T) = 0. It is not a finite H&iT) found in Ref. 
The qualitative failure of MFA might be rather unexpected, if one believes that bosons in Eq.(131) 
play the same role as phonons in the BCS superconductor. This is not the case for two reasons. 
The first one is the density sum-rule, Eq.(135), for bosons which is not applied to phonons. The 
second being that the boson self-energy is given by the divergent (at T = 0) Cooperon diagram, 
while the self-energy of phonons is finite at small coupling. 

Even at temperatures well below T c the condensed state is fundamentally different from the MFA 
ground state, because of the pairing of bosons. The latter is similar to the Cooper-like pairing of 



supracondensate 4 He atoms 1871 ] . proposed as an explanation of the small density of the single- 
particle Bose condensate in superfluid Helium-4. The pair-boson condensate should significantly 
modify the thermodynamic properties of the condensed BFM compared with the MFA predictions. 
Hence the common wisdom that at weak coupling the boson-fermion model is adequately described 
by the BCS theory is negated by our analysis beyond MFA. There is no BCS-BEC crossover in 
BFM, and the phase transition is of the BEC-type at any parameters. 



VII. CONCLUSION 

Here I have argued that attractive electron correlations, prerequisite to HTS, are caused by 
an almost unretarded electron-phonon interaction sufficient to overcome the direct Coulomb re- 
pulsion in cuprates. Low energy physics of high-temperature superconductors is that of lattice 
polarons and of real-space hole pairs dressed by phonons, i.e. mobile lattice bipolarons. Single 
polarons are thermal excitations in underdoped cuprates, but they could be degenerate in the 
boson-fermion mixture in overdoped cuprates. The superconducting state is the Bose-Einstein 3D 
condensate of bipolarons at any doping including the overdoped domain, irrespective to a pos- 
sible hybridization of single polarons and bipolarons. Our multi-polaron approach to the HTS 
problem accounts for many normal state properties of superconducting cuprates including the 



temperature-dependent spin susceptibility, nonlinear in-plane and thermally activated out-of-plane 
resistivities, the temperature-dependent Hall effect, the Nernst signal and diamagnetism above T c , 
spin and charge pseudogaps. It provides a parameter-free fit of experimental T c and describes 
isotope effects, specific heat anomalies, resistive upper critical fields, and the symmetry and space 
modulations of the order parameter. 
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